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Abstract 
Our  problem  is  to  determine  and  to  classify  by  analytic  meana  all 
planar  shock  wave  solutions  of  the  hydromagnetic  discontinuity  relations. 
The  state  ahead  of  the  shock  (i.e.,  on  the  low  density  side)  is  assumed  to 
be  known;  no  restriction  is  placed  on  the  direction  of  the  magnetic  field  in 
front.  It  is  shown,  apart  from  certain  'limit'  shocks  (e.g.,  pure  gas  shocks), 
that  the  shock  velocity  and  the  quantities  characterizing  the  state  behind 
hydromagnetic  shocks  may  be  expressed  as  simple  algebraic  functions  of  the 
discontinuity  in  the  magnetic  field  across  the  shock.  A  natural  classifica- 
tion of  all  hydromagnetic  shocks,  based  on  this  representation  of  the  state 
behind  the  shock,  is  given.  Several  useful  analytical  properties  of  the 
various  types  of  hydromagnetic  shocks  are  derived.  The  results  are  illustrated 
graphically  for  the  case  of  an  ideal  monatomic  gas.  The  relation  between 
earlier  schemes  of  classification  and  the  pi^sent  scheme  is  discussed. 
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1,     Introduction  -  Notation,  definitions,   formulation  of  the  problem 
1.1  Introduction 

In  this  report  we  shall  determine  and  classify,  by  analytic  means,  all 
'admissible  solutions'   of  the  hydromagnetic  discontinuity  relations.     We  shall 
assume  that  all  discontinuities  occur  across  a  plane  front  which  is  perpendicular 
to  the  X-axis  and  whose  motion,   if  any,   is  in  the  x-direction   [see  Figure  l] . 
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Figure  1:  Xl-cross  section  throu.:^h  planar  surface  of  discontinuity  <y  (t)  and  the 
x«axis.  i/(t)  is  parallel  to  the  (y,z)-plane}  its  motion,  if  any,  is  along  the 
X-axis.  Assuming  the  discontinuity  to  be  a  shock,  the  normal  to  «/(t)  is  always 
assumed  to  be  directed  toward  the  region  into  which  the  mass  flows.  This  region 
is  labelled  ivith  the  subscript  '1'  and  referred  to  as  the  region  'behind'  tlie 
shockj  the  remaining  half-space  is  indexed  by  'o'  and  is  referred  to  as  the  region 
in  front  of  the  shock.  The  angle  9  is  the  angle  between  n  and  H  ,  as  shown. 

Unless  otherwise  mentioned,  we  shall  assume,  in  addition,  that  the  fluid 
supporting  the  motion  is  an  ideal  polytropic  gas.  By  hydromagnetic  discon- 
tinuity relations  we  mean  the  well-known  system  of  non-relativistic  discon- 
tinuity relations  for  a  infinitely  conducting  non-dissipative  fluid  which  was 
first  derived  by  De  Hofitaann  and  Teller'-  ■'  as  a  limiting  case  from  the  corres- 

More  precise  specifications  of  the  various  expressions  appearing  in  'single 
quotes'  here  and  below  will  be  given  in  the  next  subsection. 
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ponding  relativistic  system.   Contact  discontinuity  solutions  apart,  by  an 
admissible  solution  we  mean  any  one-parameter  family  of  solutions  of  these 
discontinuity  relations  which,  for  a  fixed  value  of  the  parameter  (the  shock 
strength  parameter),  determines  the  'state  behind  the  shock',  the  'state  in 
front'  being  assumed  given. 

Let  H.  be  the  component  of  the  magnetic  field  tangent  to  the  shock 
front,  H+_  the  discontinuity  in  H.   across  the  shock  front  and  H   the 
magnitude  of  the  magnetic  field  in  front  of  the  shock.  We  shall  show  that 
when  the  ratio  h  ■  [j^+J/^o  (o^  ^^^   negative)  is  introduced  as  a  shock  strength 
parameter^  the  state  behind  compressive  magnetic  shocks  may  be  represented 
parametrically  in  terms  of  h  by  means  of  relatively  simple  algebraic  expressions » 
valid  for  all  orientations  of  the  magnetic  field  in  front.  Apart  from  its 
intrinsic  interest,  this  result  (and  others  to  be  mentioned  later)  enables  us 
to  achieve  a  more  complete  and  detailed  classification  of  the  various  shock 
wave  solutions  and  to  penetrate  further  than  has  hitherto  been  possible  into 
questions  involving  their  analytic  character. 

The  plan  of  our  report  is  as  follows:  In  the  next  subsection  we  shall 
collect  the  notation  definitions  and  equations  to  be  employed  in  the  remainder 
of  the  report  and  then  formulate  our  problem  more  precisely.  In  Section  2  we 
shall  give  a  detailed  summary  of  our  results  and  discuss  their  relationship  to 
earlier  results.  Explicitly,  we  shall  l)  exhibit  all  admissible  solutions  of 
the  basic  hydromagnetic  discontinuity  relations,  2)  classify  these  solutions  by 
analytic  means,  3)  give  several  useful  analytical  properties  of  the  solutions, 
U)  give  graphical  illustrations  of  the  solutions,  and  finally  5)  relate  our 
scheme  of  classification  and  compare  (some  of)  our  results  with  those  which  have 
appeared  in  the  literature.  The  remainder  of  the  report.  Sections  3-7}  is 
devoted  to  giving  analytic  proofs  of  the  results  collected  in  Section  2. 
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1,2  Notation,  definitional  formtLlation  of  the  problem 

Let  <^(t),  where  t  is  the  time,  be  a  surface  of  discontinuity  which 
separates  flow  regions  'u,  and  H.     of  the  fluid.  As  already  mentioned  above, 
we  shall  suppose  that  id(t)  is  a  planar  surface;  we  shall  take  it  parallel  to  the 
(y, 2)-plane  and  assume  that  its  motion,  if  any,  is  along  the  x-axis.  We  shall  assume, 
for  the  sake  of  definiteness,  that  ^  and  ^,  are  regions  of  constant  state  -  i.e.,  the 
quantities  which  desciT.be  the  flow  in  these  legions  are  independent  of  space  and 
time  variables  [see  Figure  l] .  The  discontinuity  is  therefore  uniform  in  the 
sense  that  its  description  does  not  depend  on  time. 

In  the  following  list,  all  quantities  are  assumed  to  be  expressed  in  the 
(rationalized)  MKS  Giorgi  system  of  units: 

p  s  mass  density 

'^  =  l/p:  specific  volume-i.e.,  voliime  per  unit  mass 

p  :  pressure 

e  :  specific  internal  energy-i.e.,  internal  energy  per  unit  mass 

S  !  specific  entropy  -  i.e.,  entropy  per  unit  mass 

T  :  temperature 

u  :  particle  velocity 

H  :  magnetic  intensity 

IX  :  magnetic  inductive  capacity  of  free  space 

n  :  unit  normal  to  ^ (t)  and  directed  into  region  K^ 

u  =  U'H  :  normal  flow  velocity,  or  particle  velocity  in  the  direction  of  n 

H  =  H'n  :  component  of  H  in  direction  of  n 

U  :  normal  shock  velocity  -  i.e.,  velocity  of  shock  front  in  direction  of  n 

V  =  u  -  U  :  normal  flow  velocity  relative  to  the  shock 
n   n   n  ^ 

m  =  p(u  -U  )  :  flux  of  mass  through  shock  front  -  i.e.,  mass  passing 
^     '^    through  a  unit  area  in  a  unit  time. 
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To  distinguish  between  values  of  a  quantity  Q  in    /?     and    /f,   vre  shall 
append  the  subscripts   'o'   and   '1'   and  write 

(1.1)  [qJ  =  Qi  -  Qq 

for  the  jump  in  this  quantity  across  the  shockj  here  Q  may  be  a  scalar  or  a 
vector. 

Because  of  the  conservation  of  mass  through  the  shock  front  we  have 

(1.2)  m  =  p^  V^^^  =  p^  V^^^  >  0. 

According  as  m  =  0  or  m  /^  0  we  shall  speak  of  a  shock  or  a  contact  discontinuity. 
Since  we  are  dealing  exclusively  with  one-dimensional  motion,  it  is  clear  that  n 
is  directed  along  the  positive  or  negative  x-axis.  In  the  following,  assuming 
that  m  /  0,  we  shall  employ  the  convention  that  n  points  in  the  direction  in  which 
the  fluid  crosses  ^(t).  Under  these  circumstances,  m  is  always  positive,  a 
conclusion  which  we  have  anticipated  in  (1.2).  The  region  into  which  n  points  will 
be  referred  to  as  the  region  behind  the  shock;  the  remaining  region  will  be 
called  the  region  ahead  of  or  in  front  of  the  shock.  Since  we  have  labelled  the 
region  into  which  n  is  directed  with  the  subscript  '1',  it  follows  that  the 
subscript  '1'  always  refers  to  the  region  behind  the  shock  and  the  subscript  'o' 
to  the  region  in  front  of  the  shock. 

In  terms  of  the  above  notation  the  basic  hs'-dromagnetic  shock  relations 
may  be  expressed  as  follows  jcf.  [l]-[_U],  especially  [2]  : 

^0       W  =  °' 

(^^3)  Ag  [rau  +   (p  +  nH^/2)n  -  nH^^  t]  =  0, 

A-  [m]   =0,  or  equivalently,  m[t]  -   [u^]  =  0, 

A,  m[u^/2  +  e  +  nH^  t/2]   +    [u^(p+nH^/2)-nH^T!-^  =  0, 


A^  m[s]  >  0, 


Kl^     is  the  projection  of  any  vector  Q  on  the  shock  front. 
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where  in  d^   and  ^^   the  quantities  p  and  e  have  the  form 
(1.3') 


p  =  Ap" 


e  =  p^/(y-1). 

Here,  y  is  the  so-called  adiabatic  exponent  and  i^   is  a  function  of  the  entropy 
alone. 

^   shock  wave  solution    of  A  -  A.    will  be  said  to  be  conpressive  if 
[t]  <  0,  non-conpressive  if   |r]   ■  0  and  expansive  if   [^]  >  0,      The  major  portion 
of  our  results  -  those  having  to  do  with  compressive  shocks  -  will  be  expressed 
in  the  following  notation: 

a)     ^  =  (Pi  -  Po)/Poi  *l  '  Pl/Po  =  ^  *  1» 

c)     I  =    [pI/Pq, 

e)     X  »  Y  s^Y  +  h^/2, 

f  )**  sin  9^  =  "y,  j^j»         -^°  <  »j  <  50°,      j  =  0,  1, 


(l.h) 


and 


(1.5) 


Y         1  +  s.  -/(l+s.)     -  Us.cos'^, 


£2l\      =  j 1 1 d 1 

njoy  2  cos^. 


Anticipating  somewhat,   these  solutions  are  not  necessarily  'admissible'   since 
they  may  not  satisfy  the  entropy  condition  A^. 

Since  none  of  our  formulas  will  depend  in  any  essential  way  on  the  sign  of 
cos  © .,   j  "=  0,1,  there  will  be  no  loss  of  generality  in  limiting  the  9  -  range 
to  the  interval   [-50   ,  90  ]. 
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1+s.  +  /  (1+s.)  -  Ua  .COS  9. 

2  cos  9. 
J 


Here,  H.  =  |^. |  is  the  total  magnetic  field,  and  9  is  the  angle  between  the 
J     J  "J 

positive  n-direction  and  H.  an  region   M.,  j=0,l}  ^is  the  excess  density  ratio, 

and  T,      the  excess  pressure  ratio.  The  dimensionless  ratio  h  will  be  emplojed 

as  a  measure  of  the  discontinuity  in  the  transverse  magnetic  field,  while  the 

quantity  X  will  be  used  primarily  as  a  means  of  simplifying  the  appearance  of  the 

shock  relations  when  the  present  notation  is  introduced.  The  quantity  b   .is 

n,  J 

the  Alfven  or,  as  we  shall  sometimes  say,  'intermediate'  disturbance  speed  in 
a  direction  making  an  angle  9^ =  n  H.  with  the  magnetic  field  H..   The  ratio 
s.  may  be  interpreted  as  the  squared  ratio  of  the  sound  speed 


/  -»■  v    2  -1 

to  the  Alfven  disturbance  speed  in  the  direction  of  H.,  namely  /]xH.p.     .     The 

quantities  Cg    .  and  c„    .  are,  in  Friedrichs'   terminology  ,  the  slow  and  fast 

disturbance  speeds  in    ^..     As  will  be  recalled, a.,  b      .  c      .  and  c„    .  satisfy  the 

relations     cf.    [2]   p.   12   : 

Angles  measured  clockwise  from  the  positive  n-direction  to  the  positive  Tf-direction 
are  considered  positive-otherwise,  negative   [see  Figure  Ij. 

The  existence  of  three  disturbance  speeds  was  discovered  by  Herlofson'--'  and 

\l\ 
independently  by  Van  de  Hulst'-^-'    at  about  the  same  txme.     The  speed  c     corresponds 

s 

to  the  speed  of  Van  de  Hulst's  slow  mode,  c^  to  his  fast  mode  and  b„  to  the  Alfven 

2     ^  " 

mode.  In  the  ilinit  of  small  s  =  yp/m-H  ,  Van  de  Hulst  refers  to  c  as  the  retarded 

sound  speed,  c„  as  the  modified  Alfven  wave. 
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(1.7) 


8,i        n,j    f,j 


c  .  -  a.  ^     c 


3,0   0    "f,a,     J  =0,1. 


By  the  state  P  in  front  of  a  shock  we  shall  mesn  the  set  (p  ,p  ,u  ,H  )  or  any 

__—   o  ——————————  O   O   0   o 

equivalent  set  -  e.g.,  Q  «  (p^,  s^,  sin  ©^,  u^/b^^  ^  H^).  By  the  state  f^  behind 
a  shock  we  shall  mean  the  set  |^  =  (p^,  p^  u^,  H^j  U^)  or  any  equivalent  set  -  e.g., 
r^  "  (<l^  Ip  \A>^   ^,h;m).  The  state  in  front  of  a  shock  will  always  be  assumed  known 
at  the  outset. 

Consider  now  the  basic  discontinuity  relations  (1.3)  A  ~Ai  .  Since  H 
is  continuous  across  >d(t)  and  since  f^  is  known,  A^-Ai  constitute  seven  conditions 
on  the  eight  unknowns  [T  =  ^Pi»Pi»^i>^tr  i'"^*  Clearly,  A^-A,  determine,  in 
general,  a  one-parameter  family  of  states  beliind  the  shock.  Let  e,  0  <  e  <  6^ 
denote  the  generic  shock  strength  parameter  j  e  may  be  any  scalar  quantity  -  e.g., 
W,7,h   etc.   Then  any  determination  of  the  state, 

(1.8)     \l    -  1^(6,  rj),   0  <  6<  e^, 

behind  the  shock  which  a)  depends  continuously  on  the  shock  strength  parameter  e, 
0  <  6  <  6-,  and  the  state  in  front,  and  which  b)  satisfies  the  entropy  condition 
A^  will  be  called  an  admissible  solution  of  the  discontinuity  relations. 

Our  problem  is  1)  to  determine  all  admissible  solutions  of  the  discontinuity 
relations,  2)  to  ascertain  some  of  the  general  features  of  these  solutions  and  3)  to 
classify  admissible  solutions  in  a  natural  way  by  analytic  means. 


Our  choice  of  parameters  will  be  such  that  e  =  0  corresponds  to  a  continuous 
transition  across  ^(t).  Thus,  when  we  speak  of  a  shock  we  must  require  that 
e  >  Oj  e,  may  be  infinite,  in  which  case  '<'  must  be  replaced  by  •<•  in  (1.8), 
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2,  Survey  of  resvilts 

2.1  Preliminary  remarks 

Admissible  solutions  of  the  discontinuity  relations  (1.3)  may  be 
classified  in  accordance  with  the  following  exhaustive  scheme.  First,  we 
distinguish  between  two  families  of  solutions j  contact  discontinuity  solutions 
(m«0)  and  shock  wave  solutions  (m^).  Second,  shock  wave  solutions  are  divided 
into  three  mutually  exclusive  classes s  compressive  shocks  i\2^   5  fp~  J  <0), 
non-canpressive  shocks  ( ['J^  =0)  and  expansive  shocks  ( ["ZT^  >  0).  Employing  a 
result  to  be  proved  in  Section  3>  that  only  compressive  and  non-compressive 
shocks  are  admissible  ,  we  eliminate  expansive  shocks  from  consideration.  As 
an  aside,  it  may  be  mentioned  that  this  result  is  equivalent  to  one  already  given 

r2i 

by  Friedrichs"-  ^  which  in  essence  amoujits  to   thist     A  shock  is  compressive  if  and 
only  if  the  entropy  behind  the  shock  exceeds  the  entropy  ahead   [cf .    L^]   P«   262» 
In  addition,  we  shall  show  in  Section  3  that   \j(j  <  0  Implies   [pi  >  0  and 
conversely     .     Returning  to  our  classification,  we  subdivide  the  class  of  all 
compressive  shocks  as  follows: 

Fast  magnetic  shocks  (M_  shocks):  h  sin  9     >  0;    0  <B     <  90°, 
Slow  magnetic  shocks  (M     shocks):  h  sin  8     <0;    0  <  9     <  S0° , 

S  V  o 

(?.l)  Limit  shocks   :       a)  Fast  0  -limit  shocks 

b)  Slow  0°-limit  shocks 

c)  Fast  90  -limit  shocks 

d)  Slow  90°  -  limit*** 

Here  we  have  further  restricted  the  range  of  9  to  0  <  &  <  90.  That  this 

o     -  o  ~ 

restriction  involves  no  essential  loss  of  generality,  will  be  clear  later. 


>J  ^  ** 


Expansive  shocks  violate  the  entropy  condition  (1.3 -Aj^). 
More  accurately,  [f-J  <0  implies  fpj  >  0  and  [^I  -  0  implies  [p]  ■  0, 
This  limit  is  a  contact  discontinuity  [see  subsection  2.6.^  below|. 
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The  adjective  'magnetic'  is  used  to  characterize  all  compressive  stiocks  in 

which  h  /  0.  The  use  of  'fast'  and  'slow"  will  be  justified  later.  From  the    '' 

definition  of  h  Psee  (l.Ii  -b)]  it  follows  that  H.  .  >  H    in  M_  shocks  and 

7,1    y,o     f 

H  ,<  H    in  M  shocks  .  Limit  shocks,  as  the  name  suggests,  are  shocks  which 
y,J.   y,  o     s  '  CO    , 

as  linits  of  M-C^^l  and  M  (©  )  shocks  as  0  approaches  0°  or  90°. 

Anticipating  somewhat,  this  classification  may  be  further  refined  by 
dividing  M„  and  M  shocks  into  two  subclasses  v\     ,M^   and  mI  ,  M     say, 

X        3  X     S  X       S 

according  as  the  excess  density  ratio,  excess  pressure  ratio,  etc,  are  singled-valued 

functions  of  h  for  all  allowable  values  of  h  or  are  double-valued  functions  of  h 

in  part  of  the  allowable  range.     We  shall  see  presently  how  this  division  as   to 

'type'  is  effected  analytically.     Here  we  would  like  merely  to  mention 

that  this  refinement  clarifies  the  relationship  of  W„  and  M     shocks  to  their 

f  s 

0°-limits,     As  we  shall  see  later,  Vi^^  and  M^^  shocks  have  0°-limits  which 

are  essentially  different j  the  same  applies  to  the  0°-limits  of  M^      >  M 

s     s 

shocks. 

Another  important  means  of  classifying  hydromagnetic  shocks  involves 

the  relation  of  V   ,  the  normal  flow  velocity  (in  the  region  behind  the  shock) 
n,i 

relative  to  the  shock,  to  b  , ,  the  intermediate  disturbance  speed  in  the  region 
behind  the  shock.  In  this  scheme,  a  shock  is  called  fast  if 


(2-2)    \,i  >  b^^^  ,       0  <  e  <  e^, 
intermediate  if 

and  slow  if 


> 


((.  ■     ^ — ^ — __________ 

We  have  assumed  implicitly  that  the  transverse  magnetic  field  ahead  of  and 

behind  a  magnetic  shock  is  parallel  to  the  y-axisj  we  shall  show  in  Section  k 

that  this  assumption  entails  no  loss  of  generality. 
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(2.1^)     V^^^  <  b^^^  ,        0  <  6  <  e^i 

in  (2.2)  and  (2.1|)  we  require  that  the  inequality  holds  in  at  least  some  proper 
subinterval  of  0  <  e  <  e, .  Fast  and  slow  magnetic  shocks  will  turn  out  to  be  fast 
and  slow  in  this  sense. 

Our  scheme  of  classification  is  now  essentially  complete.  That  it  is  a 
natural  one  will  become  apparent  as  we  proceed. 

In  the  following,  we  shall  exhibit  all  admissible  solutions  of  the  basic 
discontinuity  relations  and  carry  out  the  remainder  of  the  plan  detailed  at 
the  end  of  Subsection  1,2.  Contact  discontinuity  solutions,  non-compressive 
shock  solutions  (which  turn  out  to  be  'transverse'  shocks  ),  the  so-called 
parallel  shock  solutions  (characterized  by  the  condition  h  =  O)  and  perpendicular 
shock  solutions  (characterized  by  the  condition  H  =  O)  are  of  course  well  knownj 
these  will  therefore  be  accorded  a  summary  treatment. 

2,2  Contact  discontinuities 

Setting  m  =  0  in  (1.3)  we  find  that  if  H     ?^  0  then 

(2.?)  [u]   =  0,   [h]   -  0,    {>]  -  0;    [t],    [s]  arbitrary, 

while  if  H     =  0  then 
n 

(2.6)      [uj  -0,  [p  ♦  ttl^j./2j  =0i  [u^y],  [h^J,  L^],  [s]  arbitrary, 

Icf.  [2]  p.  36J,  We  conclude  from  (2.5)  and  (2.6)  that  a  shear  fLow  discon- 
tinuity is  an  inadmissible  hydrotnagnetic  motion  unless  the  magnetic  field 

normal  to  the  surface  of  discontinuity  <j(t)  vanishes.  This  interesting 

_ p— ^^ 

This  is  Friedrichs"-  •'   terminology.   De  Hofflnann  and  Tteller*-  -•   call  these 
shocks   'symmetrical'  shocks. 
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[21 
fact  was  first  observed  by  Friedrichs  '■  ■' "  and  served  as  the  basis  of  an 

investigation  by  him  and  later  by  Bazer'-  ^   of  the  motion  which  'resoires'  an 

initial  shear  flow  discontinuity  (H^  /  0), 

2.3  Non-compressive  shocks  (myfo,    [t]   °  0)-Tranaverse  shocks 

Setting   [f]    equal  to  zero  in  (1.3)  we  find  that  non-trivial  solutions 
exist  if,  and  only  if,  [cf .    [l]   -   [h]  ] 

(2.7)  m^  =  til^r'^ 

and  that  under  these  circumstances 

|uj   =0,      [p]   =  C, 

(2.8)  [uj   =  /[TF     [^J,  4^1  -  f^r^Q, 

[S]    -0. 

These  shocks  are  the  well-known  transverse  shocks,  the  word  'transverse' 
being  used  to  emphasize  that  the  jump  in  the  flow  is  along  ^(t)    .     Transverse 
shocks  may  be  termed  intermediate  in  the  sense  of  (2.3) J  for,  using  (2.7)  and 

(2.8)  we  have 

(2.9)  V     T   =  V         =  b     ,   =  b 

n,l        n,o         n,l        n,o 

Note  that  H.   ,  behind  the  shock  need  not  be  collinear  with  the  transverse 
tr,l 

magnetic  field  in  frontj  all  that  is  required  is  that  the  magnitude  of  H. 


As  seen  from  suitable  frame  of  reference,  the  flow  on  both  sides  will  be 
parallel  to  the  shock  front. 
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be  continuous  across  the  shDck.      The  magnetic  field   therefore  rotates  in  the 
plane  of  the  shock  through  some  arbitrarj"^  angle,      fsj  vanishes  because   [t] 
vanishes   [see  Section  3j • 


-vO       ooO\ 


2,1,  Fast  magnetic  shocks   (m  7^  Q.   h  sin  6^  ?►  0,  0^   /  0   ,   90  ) 

Fast  magnetic  shocks  or  M-  shocks  are  a  subclass  of  cOTipiressive  shocks. 
The  state  behind  these  shocks  is  determined  by  the  following  relations: 

|(-lh^sln9^  -  (1-^0^)  i/^^ 


M 


T,l    '^f  =  ^f 


2s^3in  9^-  (Y-l)h^ 


«f,2  ^f-  ffr*^f 


"C^h^  sin  %-(1-Sq))  t  fTThp 


2  sin  e^-  (Y-l)hj 


f  1  -  T  ^  ^f 


(2.10) 


_  l-(^/h^)sin  e^  _ 


M 


n,l 


m. 


,-1/2      ""f 


^1  n,l 


^       Po^n,0 


X 


n,l  n,l 


'?f 


[1-  Jf^Aj  sin  e^j 


TI7?* 


M 


w 


f,5  .f 


r^ 


^n,l  ^f 


^n,l 


n,l 


h.  1/2 


cos 


where 
(2.11) 


R  (h^)  =  h^ 


2       2 

Y  sin  9 


in  fe*„  I  ,2^ 

^^ 2_  -  (y-DJ.    ♦  hjsine^(2-r)(l+s^)+  ^V^^o 


-  (l-So)% 
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and  where    'f   is  used  to  identify  M„  shock  quantities. 

The  values  of  the  quantities  in  the  left-hand  sides  of  equations  (2.10) 
are  completely  determined  as  soon  as  the  state  in  front  is  prescribed  and  the 
value  of  h_  is  specified^  ►[-,  ►LA-,  Y^,  >Z^(=  rj^  +  1)  may  be  calculated  from 

"f,l  ^"^  "f,2^  ^n,l/^n,l'  ^n,o/'n,o'    Wf/^n,l  ^"^   K'^fK,!  "^^  ^«  calculated 
from  the  knowledge  of  how  7  „,   ^^/h^.  and     >?_  depend  on  h„.     These  quantities  may, 

of  course,  be  expressed  directly  in  terms  of  h„j  however,  there  is  no  real  gain 

in  doing  this. 

The  two   'types*   of  M-  shocks  -  i.e.,   the  M.       and  M^       shocks  mentioned 

above-are  defined  by  the  following  relations: 

(2.12)  m[^^     s^  >  1  -  y  sin^9^/(Y-l), 

,(2) 


(2.13)  mJ.^^     s^  <  1  -  y  sin^^/(Y-l). 


In  M^   shocks  only  the  'plus'  branches  in  equations  M„  ^  and  M„  „ 
I  i,J.      i,<f 

above  are  admissible;  the  negative  branch  leads  ultimately  to  solutions  which 
violate  the  entropy  condition.  To  be  sure,  the  plus  branch  is  not  admissible 
for  all  values  of  h-.  The  appropriate  h_  range  is: 

(2.IU)     0  <  h^  <  h-  »  -T-  sin  e^. 
I   r  Y"!     o 

In  this  range,  h„,  Y.,  V  ./b  ,,  v  /b       ,m„  and  [s]  „  vary  in  the  same  sense 
6X   X   Ti^x  xi^^       n^c  n^o  x       x 

as  h»  -  i.e.,  they  increase  or  decrease  according  as  h„  increases  or  decreases 


•»* 


Note  that  M-  , -  M„  ,  remain  unchanged  when  9     is  replaced  by  -  6  ,  and  h_  by  -h-, 

while  M„  1^  merely  changes  sign;  the  restriction  of  Q     to  the  range  0  <  9  <  90 
X ,  ^  o  o 

therefore  involves  no  essential  loss  of  generality  [see  bottom  of  p.  8]. 

In  the  following  we  shall  limit  ourselves  to  discussing  the  qualitative  behavior 
of  ^„,  Y_  m_  and  [s] „  (and  some  directly  related  quantities).  The  behavior  of  [uj 
may,  of  course,  be  obtained  from  the  knowledge  of  how  these  quantities  vary;  but, 
for  the  sake  of  brevity,  we  shall  not  discuss  the  behavior  of  this  quantity. 


-  Hi  - 


/s   _ 


As  h-  increases  from  zero  to  h_,  TLincreases  from  zero  to  2/(y-1)»  the 
familiar  gas-dynamical  limit.  At  the  same  time,  the  variables  Y^,  V^  l^n  1' 

V"^  A)   »   and   fsl-  increase  from  zero  to  infinity,  becoming  asymptotic  to  the 
n,o'  n,o'       •-  -"f 

vertical  line  of  h.  = 23in9  /(y-l).  The  dependence  of  these  quantities  on  h-  is 
clearly  illustrated  in  Figure  2  ,  where  we  have  plotted}^  ,  Y^,  etc.  versus  h^ 
for  several  values  of  the  parameter  8  >  0,  In  this  figure,  s^  =  Ij  thus  all 
curves  are  of  a  type  1,  In  Figure  3>  where  we  have  plotted  a  similar  family 
of  curves  corresponding  to  s  =  lA6,  only  those  curves  indexed  by  values  of 

0  >  37°U6'  are  of  type  1. 

(2) 
In  Klr^   shocks  both  branches  of  M_  -  and  M^  ^  ™^st  be  employed.  The 

plus  branch  yields  part  of  the  complete  admissible  solution  when  h-  is  in  the 

range 

0  <  h^  <  h^  , 

where 


^    sin  e^(2-Y)(l+s^)  +  2  cos  B^  /(Y-l)(l-8p)^+  s^Y  sin^^ 

(2.15)    h-  »  5 5 — 

^         2(y-1)  -  (Y^sin\/2) 

can  be  shown  to  exceed  h»  [see  (2.1Ii)  above].   The  remaining  part  of  the 


In  all  our  graphical  illustrations  we  have  taken  Y  to  be  5/3.  Figures  2-7 
will  be  found  at  the  end  of  this  section. 

It  is  easily  verified  that  »?_  and  Y„  remain  finite  at  h-  «  h,,  on  the  plus 

(2)  f     f  f    r 

branch  of  Ml   shocks. 
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solution  is  furnished  by  the  negative  branch  which  is  admissible  in  the 

interval 

(2.16)  ^f  -  '^f  -  ^f* 

and  which  joins  the  plus  branch  smoothly  at  h.  »  h„j  at  this  point,  R  (h-.) 
[see  (2.11)]  vanishes.  On  the  plus  branch,  Jf^,  Y^,  7^  ^/h^  ^,  V^  ^/b^  ^, 
m.  and  [s]^  vary  directly  with  h»j  as  h„  increases  from  zero  to  h„,  these 


^ 


quantities  increase  monotonically  from  zero  and  assume  finite  values  at  h„  =  h» 

at  which  point  their  h --derivatives  become  infinite.     On  the  negative  branch, 

the  quantities  i^^,  Y^,  V     -/b     ,,  V^  ^A>„  ^,m„  and   [sl-  vary  in  a  sense  opposite 
^  f       f     n,l     n,±       n,o     n,o     i  •■  -"i 

to  that  of  h_.     As  h.  decreases  f rom  1i^  to  1i„,  4„  increases  monotonically  to 

the  limiting  value  2/(y-1),  while  the  remaining  variables  increase  monotonically 

to  +00,  approaching  the  line  h     »  h„  asymptotically.     The  h„-derivatives  of  all 

quantities  are  infinite  at  h_  ■  h_j   the  joining  of  the  negative  branch  with  the 

(lower)  positive  branch  is  therefore  smooth.     These  features  are  clearly  brought 

out  in  the  type  2 -curves  of  Figure  3  -  i.e.,  those  corresponding  to  values  of 

©^  <37°U6'. 

In  all  M-  shocks,  whether  of  type  1  or  type  2  it  can  be  shown  that 

Y^,  V  tA>  -1,  V   /b   ,m-.  and  ls]»  vary  in  the  same  sense  as  Fp.  This  fact 
f '  n,l  n,l'  n,o'  n,^  f "-  •'f *■ — —  ^ 

is  illustrated  for  Y»  and  fff.  in  Figure  U  where  we  have  plotted  Y^  against  /f^. 
In  this  graph  the  upper  cxirve  corresponding  to  sin  6  =  0.1  is  a  curve  of  type 
2j  yet  Y  increases  with  increasing  ■»?„. 

We  have  already  remarked  that  M.  shocks  are  fast  in  the  sense  of  (2.2), 
The  proof  of  this  remark  follows  immediately  from  (2.10-M^  -)5  for,  since  ^^ 
and  h»  are  positive,  the  last  expression  is  clearly  greater  than  unity.  We 

therefore  have 

V^  .  >  b  T  ,    0  <  e  <  Ct  , 
n,l     n,l  *        -  1  ' 
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where  it  will  be  recalled,    e  i3  the  generic  shock  strength  parameter. 
It  can  also  be  shown  that 


(2.17)  ^n  1  '^  ^f  1  '       0  <  e  <  e^, 

and 

(2-18)  ^l^  >  Cf,o'         0  <  e  <  e^  . 

Thus,  in  M.-shocks  the  nonaal  flow  velocity  relative  to  the  shock  is  sub-fast 
behind,   and  super-fast  in  front.     This  result  is  analogous  to  the  well-known 
gas  dynamical  result;  namely,   the  normal  flow  velocity  relative  to  the  shock 
is  subsonic  behind,  and  supersonic  in  front. 
In  the  weak  shock  limit  we  find  that 

(2.19)  li-V^         =     lira/  =     c 

e->o        '  6->o      '  ' 

The  state  behind  weak  shocks  is  determined  by  the  following  set  of  equations: 

sin  0 

A\i      h^-^rrir^f  *  ••• ' 


^*^f,2         ^f  "  ^  ''f  *  •••  » 
(2.20)        AMf,3  "'f  ^1  '  Vl  '  \o  =  °f,o  *  •••  ' 


AMf,i,         Wf-     -<^f,o^f* 


r     n         ^f^f.o'^^^^o    - 

AV5     h^f' — tr^ —  Mf  *..- 

If  Jf«  is  identified  with  k(p  c,.       -  [xH  ),  where  k  is  an  arbitrary  constant, 
*x  o  1 ,o  n 

and  this  quantity  is  substituted  f or  ^„  in  (2.20),  then  the  resulting  formulas 
may  easily  be  shown  to  agree  with  fast  amaU  disturbance  fonmilas  derived  by 
Friedrichs  in  [2J    [see  p.  13J . 
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where 
(2,21) 


n,o 
2 
=f,o 


<    1. 


In  addition,  we  find  that 
(2.22)  [S]„ 


^  y(y-i)  }  ni ,  '^^o 


(l-r,)^s^ 


^V  *  oqj^). 


Here,  c  is  the  specific  heat  at  constant  volume.  Since  T„,  h_,  [u  ]»,  [%]f»  ©^c. 
depend  linearly  on  each  other  we  conclude  that  the  increase  of  entropy  across  an 
M»  shock  is  of  third  order  in  the  shock  strength. 

Finally,  we  remark  that  when  s  <  1  and  h-  =  2(l-s  )/(•,  the  values  of 
fl"-,  Y-  and  [s]»  are  independent  of  sin  9  .  Explicitly,  we  find  that 


If 


7  (^-o)  > 


(2.23) 


;p£-(l-s^)(Y-l-s^), 

0 


[s]f  =  Sr  ^°s 


2(1-s^)(y-1*s^)  +  ySq 
YS~ 


•  c^Y  log 


2(1^^)+y' 


when 

(2.2U)     hj  -  2(1-3^)/V, 

.(2) 


s  <  1. 
o 


This  property  of  Mi   shocks  is  illustrated  in  Figures  3  and  U. 


2,5  Slow  Tuagnetic  shocks  (m  /  0,  h  sin  6  <  0 .  9  /  0  ,  90  ) 

Slow  magnetic  shocks  -  that  is,  M  shocks  -  are  a  subclass  of  compressive 
shocks.  The  state  behind  these  shocks  is  determined  by  the  following  relations* 
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M     T       ^       =  h 

3^1         (S  S 


((l-3„)-rh^3in9  /2)  *  /r  (h  ) 


O       '     3 O 

(Y-l)h    +2s  sin  B^ 

3  0  O 


/ 


M 


s,2         s 


("»'' 


+  h 


((1-3^)  +  J  hgSin  e^  ±  /r  (hg) 

2  sin  e     +  (y-1)  h 
o  s 


(2.25) 


M 


n.l 


n.o 


3,3         V.3 


n,l 


"P72     b 
u  n.o 


in  -  in 

si  3 


S3    " ;p?    ^       ri+(^  A  )3inei^/2 ' 

'^1  n,l        s  o  n,o       L       »     s  oj 


M 


3>U 


L  nJi 


^n,l 


n,l 


-o 

[l^C^^A  )sinejl/2 


M 


*-  y-'s     _ 


n.l     s 


^'^       b^  T 

n,l 


„3       cos  © 
V     -  o 

n,l 


cos  &, 


[l-(Vhf)3in  9,]=^/2  , 


wliere 


(2.26)  R  (h)  -  (h  )2  ■{ 

3  3 


2   •    ^ 


-  (y-1) 


-  h  sine  (2-^•)(  1+3^) 
so  o 


2  2 

+  Us^sin  9^  +   (1-s^)    . 
0  o  o 


Here,  we  have  emplojred  the  subscript   's'   to  identify  M     shock  quantities  and 

s 

have  introduced  the  quantity 

(2.27)    \-   -^   =   ^V  "y^i^'^o 

as  a  shock  strength  parameter.  The  system  (2.25)  may  be  obtained  from  (2,10) 
on   substituting  (-h  )  for  h..  vSince  we  are  restricting  ourselves  to 
positive  values  of  9   (see  bottom  of  p.  8),  it  follows  from  the  definition 


of  M  shocks  that 
s 
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(2.28)      h  >  0. 

5 


The  two  'types'  of  M  shocks  are  defined  by  the  following  relations: 

(2.29)  M^^^    s^>l-Ysin^^, 

(2.30)  M^^^    3^  <  1  -  Y  sin^G^. 

In  M    shocks  only  the  'plus'  branches  in  equations  M  ^  and  M  „ 
above  are  admissible  and  then  only  when  h  is  in  the  range 

(2.31)  0<h  <h  =2sine. 

3  -   S  O 

In  contrast  to  the  behavior  of  the  corresponding  quantities  in  M     shocks,  the 

quantities  «"„,  T„,  V     ,A>„  t  >  ^         A>       >  m     and   [s]     do  not  depend  monotonically 
'•s'     s'     n,l     n,l'     n,o       n,o'     s  i-  -'s  '^  •' 

on  h  .     Consider,  for  example,   the  variation     of  ft   ,  J     and   fsl     with  h  .     P.s 

S  *-S    S       •-  -•  S        3 

h  increases  from  zero  to  h  ,  these  quantities  reach  finite    maxima  at 
s  s 

suitable  (in  general  not  the  same  )  interior  points  and  finally  return  to  zero 

/^ 

at  h  =  h  .  This  behavior  is  illustrated  in  Figures  5(a),  (b)  and  (e)  where 
s   s 

we  have  plotted  Y  ,  ^  and  fsi  /c  versus  h  for  several  values  of  &  >  0; 
s'  s     ^  •'s/  V        3  o  -  ' 

since  s  =1,  all  curves  are  of  type  1,  It  should  be  observed  that  the  magnitude 

of  the  transverse  magnetic  field  behind  M    shocks  is  less  than  or  equal  to 

s 

that  ahead  of  M^       shocks, 
s 

(2) 
In  M         shocks  both  the  plus  and  minus  branches  are  admissible  in 

suitable  h  -ranges.      The  plus  branch  jrields  part  of  the  complete  admissible 
s 

solutions  when  h     is  in  the  range 

It  can,   in  fact,   be  shown  that  the  value  of  h     at  which    Sl      is  maximum  exceeds 
the  value  of  h     at  which  Y     is  maximum  which  in  turn  exceeds  the  value  at  which 

__  S3 

if    is  maximum. 

'3 
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(2.32)  °  <  hg  5  hg, 


where 


^ 


-  sin  e  (2-r)(l+3   )  +  2  COS  d^  /(y-1)(1-s„)^  +  s  v^sin^^ 

(2.33)         \ 2 2 o 2 2 2_ 

3  2(y-1)  -  (y  sin\/2) 

can  be  shovm  to  exceed  h     fsie   (2.31)   above].  Tbe  remaining  part  of  the  solution 
is  furnished  by  the  minus  branch  which  is  admissible  in  the  interval 

(2.3U)         h^<    h^<h^ 


and  which  joins  the  plus  branch  at  h  =  h  j  at  this  point  ^^^^^   (see  2.26)  van- 

ishes.   Thus  IT'  <  the  state  behind,  is  a  double-valued  function  of  h  in  the  interval 
L  '  '  s 

h  <  h  <  h  .  On  the  plus  branch,  the  variables  0  ,  Y„  and  Isj  increase  to 

S  S    -      S  S         3  '-  -*  3 

finite  maxima  at  successively  larger  values  of  h     in  the  interval  0  <  h     <  h  , 

S  S    *"       3 

and  then  diminish  to  smaller  values  at  h     =  h    where  their  h„-derivatives  are 

S  3  S 

infinite.  The  negative  branch  joins  the  positive  branch  smoothly  at  h  =  h  and 

3  3 

—         —  T    1  X  ^ 

the  values  of  >?    ,  T     and    S       decrease  to  zero  as  h    decreases  from  h     to  h  . 

'S  3  *-    -"S  3  3  3 

These  features  are  illustrated  as  Figures  6(a),    (b)  and  (e)  where  we  have  plotted 

Y  .  V     and  Tsl   /c    versus  h     for  several  values  of  9^  >  0.     In  these  graphs 
s*    fs  •■  -"s    V  s  o 

s     =     lA6  and  the  type  2-curves  correspond  to  values  of  ©'q  5    30   5  the  remaining 

curves  are  type  1-curves.     In  the  range  0  <  h     <  h    the  magnitude  of  the  transverse 

magnetic  field  behind  is  at  most  the  magnitude  of  this  quantity  ahead;   in  the 

range  h    ^  h  <  h  ,    Ih     ,|    is  at  least  equal  to   1h_  ^I. 
^3        s-     s     '  y,l  y»o' 

It  should  be   clear  frcm  the  above  that,  in  contrast  to  M»  shocks,  all 

—  3       A    S 

state  variables  behind  M  shocks  are  bounded.  In  addition,  T  ,  V  ^/h     -, 

S  3    n,x   UjX 

V^  /b   ,m  and  fs]  do  not  vary  directly  with  T  as  do  T„,  V  ,/h     ,,  etc, 
n,o'  n,o'  s     '■  -'s  \s       f   n,l  n,l 

with /\ „.  In  fact,  the  graphs  of  Y  and  [s]  versus  rf  are  closed  loops.  The 
closed  full-lined  curve  (sin  Q   =  0.1)  in  Figure  h.   illustrates  this  property  for 
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the  Y  versus  7  curres, 
s       *-s 

It  is  of  interest  to  note  that  M  shocks  redace  at  h  -  h  (on  the  minus 

s  s    s 

branch  if  the  shocks  are  of  type  2)  to  what  we  shall  term  180  -intermediate 

shocks  -  that  is,  to  intermediate  shocks  in  which  the  transverse  component  of 

the  magnetic  field  behind  is  oppositely  directed  to  that  in  front.  In  the 

general  intermediate  shock,  as  we  have  seen,  the  diirection  of  the  transverse 

magnetic  field  behind  is  arbitrary. 

It  should  also  be  observed  that  the  values  of  »?  ,  I  and  fsl  are 

^s'     s  ^  -"s 

independent  of  ©     in  all  M     shocks  when  h     =  2(1-3   )/f.     In  fact,  the  values  of 
OS  so 

these  qaanti ties  at  this  value  of  h     are  given  by  the  right-hand  sides  of  (2,23). 

3 

Since  ^  and  h  are  positive  quantities,  the  denominator  in  the 

last  expression  for  V  ./h^  ^   in  (2.25)  is  greater  than  unity.  Consequently, 

n,x  n,J. 

(2-35)    V^^^  <  b;;^,  ,     0  <  e  <  e^i 

that  is,  M  shocks  are  slow  in  the  sense  of  (2.U)  above.  It  can  also  be 
s 

shown  that  .  ' 

i.e.,  M  shocks  are  super-slow  ahead,  and  that 

(2.37)      Lim  V^   -  Lijn  vj   -  c^  ^, 
6->o  "»^   e->o  ''*°         ^*° 

i.e.,  in  the  weak  shock  limit,  M  shocks  move  with  the  speed  of  slow  disturbances 

s        s 
relative  to  the  flow.  The  dependence  of  V  ,  and  V    on  h  may  be  described  as 

^  n,l     n,o     s   •' 

follows  (see  Figures  5(c),  (d)j  6(c),  (d)):   For  sufficiently  small  values 

of  h  ,  V  ,  is  sub-slow  (V  ,  <  c  ^).     fit   some  value  of  h  ,  h 
—  s'  n,l n,l    3,1  s'  a 


^8 
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say,  V  -=  c  -J  if  we  limit  oiirselves  to  the  plus  branch  we  find,  for  values  of  h 

in  the  neighborhood   of  h  but  greater  then  h  ,  that  V  ^  is  super-slow  (V  ->c  .)« 

Y    increases  directly  with  h  until  h  =  h  where  V    reaches  its  (unique)  rnaximum 
n,o  ''  3  s   s       n,o 

(1)  *   -^      (2) 

value:  In  M    shocks  V    decreases  as  h  increases  from  h  to  h  :  in  M 
s         n,o  s  ^   ^        3     s'     3 

shocks  V    decreases  as  h  increases  from  h  to  h  and  continues  to  decrease 
n.o  3  3     3 

(on  the  negative  branch)  as  h  decreases  from  h  to  h  [see  Figures  $(d)  and  6(d)3. 

S  3      3 


r  1  ^  s 

It  is  of  interest  to  note  that  S   reaches  a  maximum  at  h  .  Thus  V 

•■  -"s  s  n,o, 

m     =  p  V         and   fs]      reach  simultaneous  maxima  at  that  value  of  h    where  V^  , 
3         o  n,o  ■-  -'s 3    ■  n,l 

becomes  equal  to  the  slow  disturbance  speed  behind  the  shock.     This  property  of 
M     shocks  is  illustrated  graphically  in  Figure  7  where  we  have  plotted  V     ^/b     t, 

^'sA^n,!'    '^^-'s^'^v  against  h^  (s^  =  I/I6,  ©^  =  30°). 

There  are  two  kinds  of  weak  M     shocks,    'proper'   and   'improper.'     Vfeak 

3 

M  shocks  are  proper  if  h  is  small  when  »?  is  small  and  h  vanishes  when  ^ 
s  '^^      s  3  s  3 

vanishes.  All  other  kinds  of  weak  shocks  are  termed  improper.  The  state  behind 
proper  weak  M  shocks  may  be  obtained  from  that  behind  weak  M.  shocks  simply  by 

3  X 

replacing  r^  by  r     »  b^  ^/cf  ^,  c„       by  c_,  ^  and  h_  by(-h J  in  (2,20)  and  (2,21). 
X  s         n^os^o      x>o  s^o  X  9 

The  increase  of  entropy  across  proper  weak  M     shocks  is  of  the  third  order  in  the 

5 

shock  strength, 

/\ 
Improper  weak  M  shocks  occur  in  the  neighborhood  of  h  =»  h  ^-2   sin  9  . 
s  s    s        o 

Here, 

sin  9 

(2.38)    T_ V~  K-  ^   ^^"  %  I  *  •••'  'V  ^  ^^  %l  *"<  ^' 

1-8^-r  sin  9^ 


« 


It  is  difficTilt  to  obtain  a  simple  algebraic  expression  for  h  .  It  can  be 
shown,  however,  that  h^  is  unique  and  that  h  >  1.5  sin  9^.  The  proof  of  this 
statement,  while  not  difficult,  is  long  and  tediousj  we  shall  therefore  omit  it 
entirely. 

Graphical  analysis  indicates  that  this  is  true  for  all  allowable  values  of  h 
which  exceed  h  on  the  plus  branch  and  for  all  values  of  h  on  the  minus  branch 
in  Mg   shocks  [see  Figures  5(c)  and  6(c)]. 
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if 

(2.39)  Sq     /    1  -  Y  sin^e^  , 


and 


(2  sin  9  /y)^/^ 


0 
S  COS  5_  *~  o        s 


n  9  -  h  )^/^  +  ...»     1  »  (2  sin  9  -  h  )  >  0  , 

n  R  o  to 


(2.)^0)  -^3   =     S^TS ^2   si 

if 

(2.1a)  Sq  =  1  -  Y  sin^e^. 

In  both  cases 

(2.U2)  Yg  =  y[i  *   (Y-l)sin^9^/sJ  /j^  +    ... 

and  2 

c     Y(Y-l)sin  B 

(2J43)        [s]3  =  ^^-^ n^*  ...  . 

o 

Thus,   in  improper  weak  M     shocks  the  increase  of  entropy  across  the  shock  is 

s 

of  first  order  in  the  shock  strength  parameter  *f  . 

Finally  we  note  that  improper  weak  M  shocks  reduce  to  180  -  intermediate 
shocks  (see  p.  21)  when  Jf  is  set  equal  to  zero. 


2o6  Limit  shocks 

2.6.1  Preliminary  remarks 

Limit  shocks  are  admissible  solutions  of  the  hydromagnetic  shock 
relations  in  their  own  right.  It  is  most  instructive,  however,  to  view  them  as 
limits  of  M  (9  )  and  M  (9  )  shacks.  According  as  the  various  limit  shocks  serve 
as  lijnits  of  M„(9)  or  M  (9  )  shocks  we  have  called  them  fast  or  slow.  The 
appropriateness  of  this  nomenclature  will  become  evident  as  we  proceed. 

Note  that  the  type  relations  (2.12),  (2.13)  and  (2.29),  (2.30)  reduce, 
in  the  0  -limit, to 


(2.UU)    M^^\  M^^^  •  s^  ^  1 
m(2),  m(2)  :  s^  <  1. 
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Accordingly,  we  must  expect  two  'types'  of  0  -limits  for  M„  shocks  and  two 

'types'  of  0  -limits  for  M  shocks. 

s 

On  the  other  hand,  we  find,   employing  the  type  relations  and  the  fact 
that  3     is  intrinsically  positive,  that  the  90  -limit  shocks  are  all  of  type  1. 
Thus  we  must  expect  one  type  of  fast  90  -limit  shock  and  one  type  of  slow 
90  -limit  shock, 

2.6.2  Fast  0°-limit  shocks 

Fast  0  -limit  shocks  of  type  1  (s     >  l)t  Fast  pure  gas  shocks 
Consider  the  curves  of  Figure  2.     /s  6     approaches  zero,  the  vertical 
line  h „  =  n „  »  2  sin  ©  /(y-1),  which  bounds  all  curves  on  the  rif^ht,  moves  to 
the  left  and  coincides  in  the  limit  with  the  vertical  axis.     Clearly  then,  if 
a  non-trivial  limit  shock  exists  it  must  correspond  to  one  in  which  h„  =  0. 
It  can  be  shown  that  such  a  non-trivial  limit  shock  does  exist  and  is  in  fact 
a  gas  dynamical     shock  (abbreviated  G„  shock).     In  the  present  notation  we 
find  that  the  state  behj-nd  a  G„  shock  is  given  by  the  following: 

^f,l    ^f=  2r^f  /  [2  -  (Y-DftJ     , 


2       ^^o 

^f  ,2     '"f  =  — ? 

^1 


(2.U5) 


2-0 


2-(y-i)  n^ 


Pi 


'^f,3     Wf    "    -""^i^f' 

^f,u   tvlf  =  o> 

where  the  range  of  i^-  is 


(2.U6)  0  <7r^<  2/(y-1). 


Fast  gas  dynamical  shocks  are  fast  in  the  sense  of  (2.2)j  for,  using  G.  „  and 


the  fact  that  s  >  l.we  find  that 
o  -  ' 


i'arallel  shocks'  in  the  terminology  of  DeHoffmann  and  Teller. 
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(2.U7)     G..:  V^  >  l/^^,     0<e<e^. 

Fast  O^-limit  shocks  of  type  2  (s  <  l) :  Incomplete  smtch-on  shock  -  fast  gas 
shock  combination.  Figures  3  and  k   serve  as  aids  in  visualizing  the  limiting 

process  in  M„  shocks  of  the  present  type.  For  the  sake  of  definiteness,  consider 

Figure  3(a):  P'S   9  approaches  zero  the  upper  parts  of  the  Y»  versus  h-,-cu]rves, 

which  approach  the  line  h„  =  2  sin  9  /y-l)  asymptotically,  gradually  become 

steeper  and,  in  the  limit,  become  vertical.  The  lower  parts,  on  the  other  hand, 

approach  not  the  vertical  axis  but  a  limit-curve  -  i,e.,  the  0  curve  of 

Figure  3(a).  Thus  Y^  becomes  independent  of  h„  on  the  upper  parts  of  the  T- 

versus  h--curves  but  retains  its  dependence  on  h„  on  the  lower  parts.  Similar 

remarks  apply  to  the  «_,  V  ^/h     ,,  etc.  versus  h„-curves  of  Figure  3, 

"•I       n,±     n,l  I 

It  can  be  shown  that  the  limit  shock  is    composed    of  two  'incomplete' 
parts:  the  incomplete  switch-on  shock  (abbreviated  Sw  shock)  and  the  incomplete 
fast  gas  shock.     The  state  behind  a  G.  shock  is  determined  by  (2.U5).     Moreover, 
using  (2.U5)  G-  o>  it  is  easy  to  show  ttiat  the  gas  shock  is  fast  -  i.e., 

—  * 

if,  and  only  if,vl~  is  restricted  to  the  range 

(2,1x9)  ^(I-Sq)  ->tcrit         5    I?f     <     2/(Y-l),         s^  <  1, 

or,  alternatively,   i:^  and  only  if,  ^»  is  restricted  to  the  range 


f 


2y(1-s„) 


^2-50)    Tr^TO^  -    ^CTit      ^^f  <  "'      %<^- 


The  incomplete  G--shock  is  'completed'  by  the  incomplete  Sw  shock. 


It  is  because  of  this  limitation  on  the  e-range  that  a  fast  pure  gas  shock  is 
called  'incomplete'  when  s  <  1.  Sw  shocks  are  tenned  incomplete  for  similar 
reasons. 
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The  state  behind  an  incomplete  Sw  shock  is  determined  by: 


2 


^f  I  (1- 


^1     4  °  2'^f  i  ^^"^0^  -  ^"^"^^  *^f  /^  I  *      ^f  ^  °' 


SWg    Y^  =  Y?f' 


1  >  (Y-l)If  /2s^   , 


(,.51)     SW3     W,--b„,,jtJ/^^-//^|, 

Sw^   m^  o  pf  ^^^^  . 

The  Sw  shock,  considered  by  itself,  is  an  intermediate  shock;  for,  using 
Sw,^,  we  find: 

The  range  of  validity  of  (2.5l)  and  (2.52)  is:. 

(2.53)     0<^^<  Ti^^i^  -  2(1-s^)/(y-1),      s^<1» 

or  alternatively 

(2.5U)     0  <  y^  <  Y^^^  <  2y(1-3^)/(y-1)Sq,    3^  <  1. 

It  is  easily  verified  that  the  state  behind  the  Sw  shock  passes  continuously 
into  the  state  behind  the  incomplete  G.  shock  as  the  shock  strength  parameter 
passes  through  the  critical  value  -  e.g.,  as  »]-  passes  through 

Essentially  these  relations  were  derived  by  Friedrichs  directly  from  the  basic 
equations;  they  appeared  in  unpublished  notes  which  he  kindly  placed  at  our 
disposal.  Some  properties  of  these  solutions  were  also  mentioned  by  Friedrichs 
in  [2J.  A  more  detailed  discussion  may  be  found  in  [8j  where  it  is  shown  that  Sw 
shocks  play  an  important  role  in  the  resolution  of  an  initial  shear  flow  discontinuity. 
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'^crit  ~  '^^^-^Q^/^'^-^^'     -^^  ^^^   coraplete  range  0  <vff  -  ^A— 1»  the  ratio  Y„  increases 
monotonicaUy  vdth  ">?„.  This  is  illustrated  in  Figure  h   (see  curve  ODBCoo)  where 
we  have  plotted  Y„  against  vt„(s  =  I/I6)}   ODB  is  the  Sw  part,  BCoo  is  the  G„  part 
of  ODBCoo.  Regarded  as  functions  of  h-,  the  variables  >f»,  Y„,  m  ,  [s]  „,  on  the  Sw 

branch,  are  double-valued,   [see  Figure  3»] 

f       f 
The  relationships  of  V  ,  and  V    behind  and  ahead  of  the  switch-on 

n,i     n,o 

shock— fast  gas  shock  (abbreviated  SwG.)  combination  to  the  small  disturbance 
speeds  ahead  of  and  behind  the  shock  are  given  by  the  following: 


(2.5$)    SwG^: 


^n,l  "  ^n,l  <  '=f,l'   ^^V'^-^t. 


(2.56)    SwG^:   V^^^  >  b^^^  >  a^,   0  <  e  <  e^. 

The  state  behind  weak  switch-on  shocks  is  determined  by: 

ASw^    i^  =   2(1-3^)  71^-*-   ,.., 
A  Sw^    Y^  =  Y  »[/•..  > 
(2.57)   ASW3   Wf=  -t>^n,o  V  •••' 

aSw^    vf  t  =  V^  „  =  b^  ^  +  ...  . 
**  5     n,l    n,o    n,o 

In  addition,  we  find: 

r,    %  f  (t-1)(1-3.)  _2         0  Y(Y-l)hj! 

(2.58)       Hf  =  -^^-5r -"f*--    8,(1-3    

O  0     0 

Thus,  the  rise  in  entropy  across  an  Sw  shock  is  of  the  second  order  in  the  excess 
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density  or  pressure  ratio  and  of  the  fourth  order  in  the  variable  h_. 

2.6.3  Slow  0°-lijnit  shocks 

Slow  0°-limit  of  type  1  (s     >  1):  Continuous  transition 

In  the  9     =0°  limit  all  discontinuities  of  the  M^       shock  vanish, 
o  s 

producing  a  continuous  transition  across  the  surface    o(t).     The  liiniting 

process  is  most  easily  visualized  with  the  aid  of  curves  of  Figure  5.     In 

Figure  5(a),   for  example,   the  ranges  of  Y     and  of  h     decrease  with  decreasing  9 

s         s  ^ 

and  reduce  in  the  0  -limit  to  zero. 

Slow  0  -limit  shocks  of  type  2  (s  <  1):  Slow  gas  shock— -switch-on  shock 
combination. 

As  was  the  case  with  Mi.   shocks,  slow  0  -limit  shocks  of  type  2 

consist  of  two  incomplete  parts.  One  part  is  an  incomplete  slow  gas  shock 

(abbreviated  G  shock);  the  other  part  is  an  incomplete  switch-on  shock. 

The  state  behind  the  incomplete  G  shock  is  determined  by  (2.h5)  with 

s 

'f   replaced  by  's'.     It  is  slow  because  ^     is  restricted  to  the  range 

(2.59)  0  <73  5Jfcrit     =     2(1-s^)/(y-1), 

[cf.  (2.53)]  and  because  in  this  range  (and  in  this  range  only) 

(2.60)  v^^^  <  b^^,  , 

the  equality  holding  only  when/^   =  '^crit* 

The  state  behind  the  Sw  shock  is  determined  by  (2.5l)-(2.5U)  with  h^ 
replaced  by  h  .  This  Sw  shock  differs  from  the  one  in  the  SwG^  combination 
in  one  respect  onlyj  here,  H  ^  is  negative  vrhereas  H  ^  exceeds  zero  in  the 
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SwG«  combination. 

In  the  (*f  ,  Y  ) -plane  the  graph  of  Y  against  ^  is  a  closed  loop 

S    S  So 

[see  Figure  U] .  This  loop  consists  of  the  arc  OAB  (G  shock)  and  the  arc  BDO 

*•  3 

(2) 
(Sw  shock).  The  smooth,  closed,  nearby  loop  is  an  M    curve  corresponding  to 

s 

a  small  value  of  9  ,  namely,  ©  -  arcsin(0,l). 

The  limiting  process  is  graphically  portrayed  in  Figure  6.  Consider 

Figure  6(a)  for  example!  These  curves  may  be  divided  into  two  arcs,  one  in 

the  range  0  <  h  <  h^(e^),  the  other  in  the  range  h(&^)  <  h  <  h  ,  where  h  (©  ) 
sso  os"~s       so 

is  the  value  of  h  for  which  ?  is  maximum,  say.  As  6  approaches  zero  the 
s  s  s 

first  arc  becomes  steeper  and  in  the  limit  coincides  with  part  of  the  Y  -axis. 

The  remaining  part  reduces  to  the  switch-on  shock. 

The  relationship  of  V^  ,  and  V^   behind  and  ahead  of  the  slow  gas 
^     n,l     n,o  ^ 

shock-switch-on  shock  combination  (abbreviated  G  Sw  combination)  is  given  by 
the  following: 

(2.61)  GSw:  V°  T  <  min(a?,  b^  t),   0  <  ^i  <  ^  .^    (on  G  branch), 

8     n,l  -     In,!'       '   '^crit        s 

^,1  ■  ^n,l»   °  ^  ^g  ^  ^crLt  ^°"  ^^  ^^^"''h^' 

(2.62)  QgSwt  ^n  o  '^  ^'o*      ^  ^  ^s  -  ^cidt  *  ^°"  "^s  ^"^  ^"  branches)  , 

g 

V    is  at  most  b    when 
n,o  n.o 


The  basic  shock  relations  do  not  uniquely  determine  the  direction  of  the  trans- 
verse magnetic  field.  If,  however,  Sw  shocks  are  regarded  as  0°-limits  of 
M-(9  )  and  M  (6  )  shocks,  then  it  follows  (since  we  are  always  assuming  9  >  O) 
that  H  T  >  0  in  fast  Sw  shocks  and  H  ^  <  0  in  slow  Sw  shocks. 
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(2.63)  0  <  "^o  <    2(l-s   )/(2s     +  Y  -1)  ,         (on  G     branch). 


and  at  least  b         in  the  ranges 
n>o 

2(l-s^)/(2s^-H  Y-1)   <7s^7crit     '  (on  G^  branch)      , 

(2.6h) 

0<7<'?      •+>  (onSw  branch) 

2,6,U  Fast  90°-limit  shocks  (perpendicular  shocks) 

Fast  90  -limit  shocks  are,    in  the  terminology  of  de  Hoffmann  and 
Teller,  the  so-called  perpendicular  shocks.     We  have  seen  that  90  -limit 
shocks  are  of  type  1.     The  state  behind  fast,  perpendicular  shocks  is  determined 
by  the  following  relations: 


'f,l      ^f  "  ^f    » 

_        Y  ^^[l  +  (Y-1)  ^f  AsJ 
[l  -  (y-1)  ^^/2\ 


Ko     ^r'-^ -r^  »       o<?<2/(y.i), 


^^•^^^  f  r  ,       s  -   /  il/2 


p.  .      V 


f  <o       ^      [1^3^^  (2^)  V2]^/^  ,     ./?7r 


[i+Sq+(2-y)  17^/2] 


1/2 


It  follows  from  P^  ,  that  v    varies  in  the  same  sense  as  w...  From 
f,3      n,o  1 


this  we  conclude  that 


(2.66)     V^   >  b„  S^^^f^'^     -  (^n  o  ♦  ^o^^^^  =  °f  o- 
n,o    n,o    0         n,o    o         i,o 
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Comparing  the  expression  for  c-  ,   -  i.e., 

1/2 


b^  ^     j  i-»-3^[h-(y-^i)  7i^/2j  *  n^{3--f)/2  *  ^/(2wy)[i-(y-i)7?^A] 

f  ,1  "  "^    1  1  .  (^.1)  ;^^/2 


with  the  expression  for  V  ^  in  P„  -  we  find,  in  addition 

n,j.    i,^ 

(2.67)  V^^i  <  c^^i  ,  0  <  e  <  c^. 

Weak  fast  perpendicular  shock  equations  may  be  derived  from  (2,20)- 
(2.22)  simply  by  setting  ©^  "  90  and  r^  -  b^  Jz^  ^  -  0  in  these  equations. 

2.6.5  Slow  90^-limit  (contact  discontinuity) 

The  w'^-limit  of  an  M^"^^  shocks  is  a  special  case  of  the  contact 

s 

discontinuity  solution  characterized  in  (2.6).  /\s  in  (2,6)  we  have 

(2.68)  m^  =0,  [ujg=0,  [p  +  nH^y2]^=  0,   O^^J^ arbitrary. 

In  contrast  to  (2.6),  however,  M  and  [s]  are  no  longer  arbitrary.  For 

s      s 

example,  it  can  be  shown  that 
(2-h  )h 

(2-^9^         %'   ^s\  (r-l)h    • 

O  S 

Since   [s]    /c     =  log(Y  /t[^),    [s]     also  depends  on  h  . 

S   V  S   S       S  S         . 

2,7  Concluding  remarics  -  Survey  of  the  literature 
2,7.1  Concluding  remarks 

We  add  here  several  remarks  with  the  object  of  completing  the 
qualitative  picture  of  hydromagnetic  shocks.  The  first  is  that  the  excess 
density  ratio  in  all  compressive  shocks  never  exceeds  the  gas  dynamical  upper 
limit  2/(r-l)«  This  upper  limit  is  actually  assumed  in  M^  shocks;  in  M^  shocks 
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the  maximum  value  of  the  excess  density  ratio  is  less  than  2/(y-1).     The 
second  remark  may  be  expressed  as  follows:     The  pressure  behind  a  magnetic 
shock  corresponding  to  a  given  value  of  the  excess  density  ratio  isj  owing  to 
the  presence  of  transverse  magnetic  fields,   greater  than  the  pressure  behind  a 
pure  gas  shock  at  that  value  of  the  excess  density  ratio    [see  Figure  iTj .      This 
result  is  an  immediate  consequence  of  the   relation    [cf.   Lust    [hjaj 

(2.70)         Y  =  -i]L3 .    rir-l)h^^   _  ^^^^  2/(y-i)  , 

2-(y-1)  71  2s^[2-(y-1)  n  ] 

which  is  the  hydromagnetic  analog  of  the  gas  dynamical  Hugoniot  relation.  The 
third  remark  involves  a  result  already  implicit  in  the  above  development  but  which 
nevertheless  deserves  to  be  singled  out  for  special  emphasis.  It  is:  the  magnetic 
field  behind  compressive  shocks  is  finite  for  all  values  of  the  shock  strength 
parameter.  Finally  we  remark  that  the  entropy  behind  a  compressive  shock 
varies  in  the  same  sense  as  the  mass  flux  or,  alternatively,  as  the  normal  flow 
velocity  (relative  to  the  shock)  in  the  region  ahead  of  the  shock.  This  result 
follows  directly  from  the  differential  relation  [see  Subsection  3,2} 


(2.71)   T^d[s]  =  I  [q2  .  [rn  f/H^ 


mdm, 


which  we  believe  has  not  appeared  in  the  literature  thus  far.  It  is  valid  for 
all  compressive  shocks;  the  last  term  is  absent  in  pure  gas  shocks  and  in 
perpendicular  shocks. 

2.7.2  Survey  of  the  literature 

The  study  of  compressive  hydromagnetic  shocks  of  finite  amplitude  was 
initiated  by  de  Hoffmann  and  Teller''  J,  who  as  already  mentioned,  derived  the  hydro- 
magnetic  analogs  (relativistic  and  non-relativistic)  of  the  Rankine-Hugoniot 
relations  of  gas  dynamics.  Shock  wave  solutions  of  these  relations,  corresponding 
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to  special  orientations  of  the  magnetic  fields  in  front  of  and  behind  the 
shock  front-  explicitly,   parallel,  perpendicular  and  syinmetrical  shocks 
(i.e.,  our  intermediate  shocks)- were  discussed  in  some  detail  by  these  authors 
and  the  extreme  i^lativistic  and  non-relatiristic  behavior  was  examined.      The 

steady  shock  wave  solutions  of  the  non-relativistic  set  of  hjd  romagnetic 

[31 
discontinuity  relations  was  then  taken  up  by  Heifer '■■'   who  discussed  shock  wave 

solutions  corresponding  to  more  general  orientations  of  the  magnetic  field 

(i.e.,  oblique  shocks).     Heifer's  analysis  involves,   in  addition  to  T  -  p-./p  , 

the  three  parameters  Q,   s  and  x,  which  in  our  notation  and  system  of  units 

[see  Subsection  1.2]   may  be  expressed  as  follows: 


«-^ 


(2.72)  s  ■=  tan  9     ,  0  <  ©     <  90°  , 

o  o  ' 


X  =  tan  9, /tan  9     . 
1  o 


Heifer  divides  magnetic  shocks  into  three  classes     a)x>l,  b)0<x<l, 
c)  X  <  0,  and  derives  bounds  on  s  and  x  which  serve  to  delimit  the  ranges  of 
s  and  X  in  which  the  shocks  are  admissible.     Heifer's  results  for  the  case 
Y  =  5/3  are  summarized  graphically  in  Figures  3-7  of  his  paper  where  he  plots 
Q  =  Q(Y,  X,   s)  versus  Y  for  fixed  s  and  for  several  values  of  x.      Using  Heifer's 
results  one  could  obtain  a  plot  of  Y  vers\:is  x  for  several  values  of  s  and  fixed 
Q.     It  would  then  perhaps  be  possible  to  compare  our  results  with  his.     Such  a 
procedure  would,   however,  take  us  too  far  afield.     We  shall,  instead,  be  satisfied 
with  noting  the  following  relations  betvreen  Heifer's  magnetic  shocks  and  ours: 
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a)  X  >  1 


(2.73) 


b)  0<5C<1 


M 


M 


M 


M 


(1) 

f     ' 

(2) 
f     ' 

(1) 

s     ' 

(2) 


0  <  hj.  <  h^ 


2  sin  9„/(y-1), 


^ 


0  <  hj  <  h^,    [cf.    (2.15)] J 


0  <  h     <  sin  ©^  , 
s  -  o  ' 


0  <  h     <  sin  ©     J 
s  ~  o  ' 


c)  X  <  0 


M^^,^     sin  ©     <  h     <f  ,    [cf.    (2.33)]    . 

S-"  OSS'-  -■ 

[U(b)]  ,, 


The  work  of  Hist •-        -'    is  mainly  numerical  in  character.     LUst's 

analysis  is  based  on  the  fact  that  the  magnetic  shock  velocity  and  all  quantities 

behind  may  be  calculated,  for  a  given  state  in  front,  once  a  suitably  defined 

Mach  number  is  prescribed  and  once  the  dependence  of  the  density  ratio  on 

this  Mach  niamber  is  known.  Ltlst  shows  that  the  density  ratio  satisfies  a 

cubic  equation  whose  coefficients  are  cubic  polynomials  in  the  square  of  the 

Mach  number.  He  then  solves  this  equation  numerically  for  the  case  of  y  =  5/3 

and  determines  the  velocity  of  the  shock  and  the  state  behind  the  shock.  Because 

our  choices  of  shock  strength  parameters  differ  from  those  Liist  has  employed 

(i.e.  U  /c„,  U  /c  ,  U  A>  ))  a  direct  comparison  of  numerical  results  is  again 

not  possible.  Nevertheless,  by  making  use  of  the  relations  V    >  c»  „ 

'  "^  n,o    f,o 

(in  M„  shocks),  V    =  b    (in  intermediate  shocks),  v    >  c    (in  M 


n.o 


n,o 


n.o 


s,o 


shocks)    [cf.    (2.18),    (2.9),    (2.36)]   it  is  possible  to  identify  our 

M  ,  /IfVen  and  M     shocks  with  the  three  tjrpes  of  shocks  distinguished  by  LUst 


I  see    [U(b)]   p,  13Ul,     namely,  those  in  which  lijii  V 

e->o   ' 


=  c„  ,  those  in  which 
o    f  ,o' 


n,o 


=  b 


n,o 


and  those  in  which  lim  V^ 
e->o 


n.o 


=  c. 


s.o 


This  identification  having 


been  made,  it  appears  that  there  is  complete  agreement  on  a  qualitative  level 
between  LQst's  numerical  results  and  our  analytical  and  numerical  results. 

The  inspiration  for  the  present  work  comes  from  the  report  of  Friedrichs' 


L2] 
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Friedrichs   (as  does  Liist)   finds  three  types  of  shocks,   the  so-called  fast  , 
intermediate  and  slow,  which,   in  the  weak  shock  limit,  move  with  the  fast, 
interrnediate  and  slow  disturbance  speeds  respectively.     In  contrast  to  LUst's 
investigation,  however,    Friedrichs'  work  enjojrs  the  advantage  of  being  chiefly- 
analytical  in  character.     Our  M„  shocks  correspond  to  Friedrichs'    fast  shocks. 
This  is  to  be  expected  since  the  essential  features  of  his  fast  shocks, 

namelv  1h     ,|>1h       l.V     ,>b     ,,   have  been  built  into  our  definition  of 
iioiUDj-jr    I   y^Q^i        I   y,o''     n,l         n,l' 

M     shocks.     Our  M     shocks  are  equivalent  to  Friedrichs'    slow  shocks  provided  h 

f  3  ° 

is  restricted  to  the  range, 

C  <  h  <  sin  ©„  . 
s  -     o 

It  is  important  to  note,  however,  that  Friedrichs  limits  his  slow  shock  entirely 
to  this  range  and  finds  as  a  consequence, 

V^  T  <  c  T 
n,l  -   s,l 

This  limitation  of  the  range  of  h  is  unnecessarily  restrictive j  for,  as  we  have 

(1)  ^ 

seen,  the  appropriate  range  is  0  <  h  <  2  sin  9  in  M    shocks  and  0  <  h  <  h 

S  ^        OS  s     s 

(2)  ^  r  n 

in  ir   shocks  where  h  >  2  sin  9       see  Subsection  2,51.  While  it  is  true  that 
s  sol-  -^ 

V^  ^  <  c  T  in  the  range  0  <  h  <  sin  9  ,  it  is  not  true  throughout  the  complete 

admissible  range  at  h  .  Incidentally,  this  resolves  a  discrepancy  between  a 

result  of  Lttst  and  one  due  to  Friedrichs,  LUst  shows,  by  numerical  means,  that 

slow  waves  exist  in  the  region  b  t^V  t>c  ,  and,  more  important  still,  that 

the  slow  wave  is  connected  to  the  180  -intermediate  shock  [see  Subsection  2,5] . 

The  fact  that  180  -intermediate  shocks  are  not  so  connected  in  Friedrichs' 

analysis  is  due  to  his  unnecessary  limitation  of  the  h-  range. 

s 

Our  analysis  differs  from  Friedrichs'  in  several  other  respects.  First, 
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it  should  be  observed  that  J-riedrichs  specifies  the  average  state  (that  is,  the 
average  of  like  quantities  on  both  sides  of  the  shock) j  the  shock  velocity  and 
the  state  behind  are  then  determined  in  terms  of  the  average  state  and  a  suit- 
able shock  strength  parameter.  This  type  of  solution  of  the  shock  relations 
has  the  advantage  of  simplifying  the  analysis  and  of  leading  directly  to  many 
important  properties  of  the  various  shocks.  In  initial  value  problems,  howe-^er, 
it  i?  the  state  in  front  that  is  usually  knownj  shock  v/ave  solutions  in  which  the 
state  in  front  is  specified  are  therefore  preferable  in  this  corjiection. 

Basically,  however,  the  advantage  of  our  procedure  stems  from  the  fact 
that,  given  the  state  in  front,  we  have  been  able  to  determine  the  state  behind 
magnetic  shocks  in  terms  of  relatively  simple  algebraic  functions  of  the  jump  in 
the  magnetic  field.  By  varying  the  state  in  front  we  are  thus  able  to  obtain  a 
comparatively  complete  picture  of  the  family  of  all  shock  wave  solution  of  the 
basic  discontinuity  relations.  In  addition,  important  properties  of  these  solutions 
may  readily  be  derived.  Thus,  for  example,  as  we  have  seen,  the  magnetic  field 
behind  magnetic  shocks  is  finite  for  all  values  of  the  shock  strength  parameter? 
an  analytic  proof  of  this  fact  seems  never  to  have  been  given  before.  Our 
particular  representation  also  enables  us  to  obtain  a  more  complete  classification 
of  the  various  shock  wave  solutions.  For  instance,  our  classification  of  magnetic 
shocks  according  to  'type'  (see  Sections  2.U  and  2,5)  may  be  regarded  as  a 
refinement  of  the  scheme  of  classification  given  by  Friedrichs,  As  we  have 
seen,  this  refinement  has  the  virtue  of  clarifying  the  connection  of  the  slow 
and  fast  magnetic  shocks  with  their  zero-degree  limits. 

Finally,  it  should  be  mentioned  that  our  definitions  of  fast  and  slow 
shocks  differ  somewhat  from  the  corresponding  definitions  given  by  Friedrichs. 
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We  define  a  shock  to  be  fast  ifu  ,-U  >b  ,,  slow  ifu  ,-U  <b  ,, 

n,l    n    n,l'         n,l    n  -  n,l' 

and,  in  addition,  we  require  the  inequality  to  hold  in  at  least  a  proper 
subrange  of  the  shock  strength  parameter  range;  Friedrichs  does  not  include 
this  last  requirement.  Our  definitions  allow  us   to  regard  the  0  -limits  (see 
Sections  2.6.2  and  2.6.3)  of  slow  and  fast  shocks  as  slow  and  fast,  respectively. 
These  definitions  and  the  'type'  relations  enable  us  to  place  0  -limit  shocks 
in  their  proper  context  within  the  family  of  solutions  obtair«d  by  varjdng  the 
direction  of  the  magnetic  field  in  front. 
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Figure  2 

Post    mognetic  shock   C  s^^  .  I ) 

Al I  curves   ore  type  I 
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Figure  2,  Illustrating  the  dependence  in  M«  shocks  of  Y^,  *2.«,  7  ^/b  ^ 

°x.  nA>  i>  ^  „A>  ^j  [s]^/-  on  the  shock  strength  parameter  h-  for 
x«^  Ti^^       n^o  n^o     x/c  x 

several  values  of  the  parameter  9  .   Here  s„  =  1  and  9  varies  from 

0  0  0 

He  -  90  .  ^11  curves  are  of  type  1.  The  maximiim  value  of  h^  is 
h »  =  2  sin  ©  /y-l  =  3  sin  9  .  y\s  h_  increases  to  h»,  t^  increases 
monotonically  to  2/y-1  =  3j  the  remaining  dependent  variables  increase 
monotonically  to  infinity,  approaching  the  line  h^  =  h„  asymptotically. 
As  ©  approaches  zero  the  asymptotes  and  their  associated  curves  approach 
the  vertical  axis.  The  0  -limit  curves  are  independent  of  h-  and  corres- 
pond to  the  fast  pure  gas  shock  (see  Section  2.6,2), 
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Figure    3 


Fast  magnetic  shock  (Sq  =  — ) 


Type  I:    9^>  37"46' 
Type  2:  8„<  37*  46' 


3sin90° 
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Figure  3'     Illustrating  the  dependence  in  M„  shocks  of  Y  ,  /{   ,  V     ^/b     ,, 

0-  ^/b     T >  V       /b       ,    [s]^/      on  the  shock  strength  parameter  h„  for  several 
if-i-     ^j-*-       n^o     n^o  L/C  I 

values  of  ft   .     Here,   s     =  lA6j   curves  corresponding  to  values  of  9     <  37  h6' 
are  type  2  curvesj  9     >  37°U6',   type  1  curves.     The  C   -curves  represent  the 
state  behind  an  incomplete  switch-on  shock.     The  maximum  value  of  h„  in  M). 
shocks  is  hj  =  2  sin  9  /(y-l)   =  3  sin  9  :  in  M»       shocks,  the  maximum  value 
is  h-  >  h-  [see  (2.1^)].     As  h-  approaches  h-  and  ^-  approaches  2/(y-1)  =  3> 
the  remaining  dependent  variable  becomes  infinite,   approaching  h-  =  h» 
asymptotically,     /s  9     approaches  zero,   the  upper  parts  of  all  curves  become 
steeper  and  in  the  limit  become  vertical.     The  lower  parts,  on  the  other 
hand,   approach  not  the  vertical  axis  but  the  incomplete  switch-on  shock. 
The  limiting  solution  is  the  incomplete  gas  shock— switch-on  shock 
combination  (see  Section  2.6.2  and  Figure  U). 
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Figure      4 
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Figure  U»  Illustrating  the  relationship  between  Y  and  "^  in  SviG  and  G^Sw 

(9)  (?) 

shocks  and  in  neighboring  M),  '^  and  M^^  shocks.  The  curve  0;^BCoo  is  a 

X         s 

pure  gas  shock;  the  arc  0/B  is  the  incomplete  G  shock;  the  remaining  arc 
BCoo  is  the  incomplete  G-  shock.  The  arc  ODB  represents  the  incomplete 

switch-on  shockj  the  closed  loop  O/.BDO  represents  the  G  Sw  combination; 

(2)  (2) 

the  curve  ODBCoo  the  SwG„  combination.  The  M^  '^  and  M^  '  shocks  reduce 

in  the  0°-limit  to  G  Sw  and  SwG„  shocks,  respectively. 
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Figure    5 
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Figiire  5»     Illustrating  the  dependence  in  M     shocks  of  Y  ,  ■>f  ,  7     -tA^„  t> 

c     ,  A>^  -■  J  V^     A>       »    fsl    /c     on  the  shock  strength  parameter  h     for  several 
s,l     n,l'     n,o'    n,o'    ^  -^s     r  b       ^  g 

values  of  the  parameter  ©  .     Here,  s     =1;  all  curves  are  of  type  1,     The 

(l)  ^  ^ 

maximuir,  value  of  h     in  M^       shocks  is  h     =2  sin  9^.     As  h     approaches  h^ 
as  s  OS  s 

the  M         shock  reduces  to  a  180   -intermediate  shock  -  i.e.,  to  an 
s 

intermediate  shock  in  which  the  transverse  magnetic  field  behind  is  directed 

oppositely  to  that  ahead.     In  the  0°  limit,   all  discontinuities  of  the  M^ 

s 

shock  vanish. 
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Figure   6 


Slow  magnetic  shock      (Sn=^) 


0°   shocks 
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Figure  6«  Illustrating  the  dependence  in  M  shocks  ofY,»,V  ^/b  -, 

c  ,A>^  n»  V   /b   ,  fsl  /c  on  the  shock  strength  parameter  h  for 
3,1'  n,l'  n,o'  n,o'  >-  -'s  v  '^   ^        s 

several  values  of  the  parameter  9  .  The  curves  corresponding  to  values 

of  9  <  30  are  type  2  curves j  60  -curves  are  type  1  curves.  The  0  -curves 

represent  the  state  behind  an  incomplete  switch-on  shock.  In  type  1  curves, 

the  maximum  value  of  h  is  h  =2  sin  9  :  in  type  2  curves  the  maximum 

s  s  o'  *^^ 

value  of  h     is  h     [see   (2.33)].     In  the  0  -limit  the  curves  reduce  to  those 
corresponding  to  the  slow  gas  shock— switch-on  shock  combination  [see  Section 
2,6.3,   also    Figure  h\.     Those     portions  of  the  vertical  axis  intercepted 
by  the  switch-on  shock  curves  correspond  to  the  slow  pure  gas  shock. 
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Figure  7 
Slow  magnetic  shock 

t^o=r5.«o  =  ^^'^ 

[Sjg/Cy    vs  h. 


.200- 


s- 


-19  ' 


Fisure__7.  Depicting  the  variation  of  V^  ^/h^  ^,   V^  ^/b^  ^,   Cg^i^n,!  ^^ 

[s]  /c  with  the  shock  strength  parameter  h  ,  Note  that  [s]g/c^  and  V^  ^Aj^  ^ 

s  •  ^ 

reach  their  maximum  values  at  the  point  where  V  ^=  c  .,  -  i.e.,  at  h  =  h  . 

n,  J.       3  9  *^  IS 

The  plus  and  minus  signs  which  appear  as  superscripts  (e.g.,  ^'^f  ^n'l*  ^^^*' 
denote  the  plus  and  minus  branches  of  the  curves.     Note  that  V       /b         is 
unity  when  h     •=  1i     =  2  sin  B     =  1.0   and  when  h     =  sin  9     -  0,5. 


-  so  - 


3,  The  inadmissibility  of  expansive  shocks  -  The  relation  between  d[s]  and  dm 
in  compressive  shocks 

3.1  The  inadmissibility  of  expansive  shocks  and  related  results 

In  this  subsection  we  shall  show  that  only  non-expansive  shocks  satisfy 
the  basic  discontinuity  relations  A,  -  A^  of  (1.3).  Our  method  of  proof  will 
consist  in  showing  that  i)  [f]  <  0  implies  [s]  >  0,  ii)  \T]   =0  implies 
[3}  =  0  (and  conversely)  and  iii)  [^J  >  0  implies  [s]  <  0.  As  a  by-product 
of  our  proof  of  i)-iii)  we  shall  find  that  [t]  <  0  implies  [p]  >  0  and  con- 
versely, and  that  [^1  vanishes  if,  and  only  if,  [p]  vanishes.  Note  that  expansive 
shocks  violate  (1,3  -A^) because  of  iii)  and  hence  are  inadmissible.  Note  also 
that  statements  i)  -  iii)  are  equivalent  to  the  following  assertions: 

a)  The  entropy  behind  a  shock  exceeds  the  entropy  in  front  if,  and  only 
if,  the  density  behind  exceeds  the  density  in  front j  b)  the  entropy  behind  a 
shock  is  the  same  as  that  in  front  if,  and  only  if,  the  density  behind  is  the  same 
as  that  in  front. 

In  proving  these  results  we  shall  require  only  that  the  caloric  equation 
of  state  of  the  medium,  namely, 

(3.1)      P  =  g(^,S), 

satisfies  the  following  conditions: 

(3.«   c^x  <  0.  ( J)^  >  0,  (i)^  >  0 , 

and  that 


C3.3)    (11)^   <  0, 


here  p  is  to  be  regarded  as  a  function  of  f  and  e. 


1^1   denotes  the  partial  derivative  of  x[^   X(Y,Z)]  with  respect  to  Y,  Z 


being  held  fixed. 
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For  any  value  of  S,  the  first  two  relations  of  (3.2)  require  that  the  function 
gCt^S)  is  convex  downward.  The  last  relation  of  (3.2)  requires  that,  for 
constant  specific  volume,  the  pressure  increases  with  entrDpy.  In  (3.3)  it  is 
assumed  that,  for  constant  specific  internal  energy,  the  pressure  decreases  as 
the  specific  volume  increases,  ^ny  polytropic  ideal  gas  satisfies  (3.2)  and 
(3.3)  since 


(3.a) 


(11)3"-?.      (f^]-^'-^'# 


where  the  constant  c  is  the  specific  heat  at  constant  volume. 

^e  begin  our  proof  of  statement  i)  -  that  is,  [t]  <  0  implies 
[s]  >  0,  by  introducing  the  hydromagnetic  analog  of  the  Hugoniot  relation: 


(3.5)      [e]  '  -  W 


?*^      I    ,  P-(P,  *P„)/2 


This  relation,  which  was  first  derived  by  Ltist '•''•'    ,  may  be  obtained  from  the 
energy  relation  Ai  after  eliminating  the  velocity  terms  by  means  of  A^  -  A_ 
and  cancelling  the  factor  m  ^  0. 

Let  (  T  ,  p  )  and  (  f, ,  p^ )  be  the  specific  volume  and  pressure  in  front 
of  and  behind  the  shock  respectively,  e  and  e,  the  corresponding  specific 
internal  energies  and 


(3.6) 

%' 

p   -    h(r  ,6^) 

and 

(3.7) 

h'- 

p     =     h(r,e^) 

the  curves  of  constant  specific  internal  energy  through  ("^qsPq)  and  C^njPn). 
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Our  first  step  will  be  to  prove  that  K. ,  in  a  compressive  shock,  lies  above  E 
and  that  p^  exceeds  p  .  To  this  end  we  note,  employing  (3.5),  that 


(3.8)  ®l  -  ®o  ^  ° 

in  a  compressive  shock.  Next,  we  observe  that 

(3.9)  (II)  >  0. 

This  relation  follows  from  the  identity 

after  setting  d'^  =  0  and  making  use  of  the  last  assumption  of  (3.2).  From 
(3.8)  and  (3,9)  it  follows  that  K,  lies  above  E  in  the  ( f ,p)-plane.  This 
result  and  (3.3)  then  imply  that  p,  exceeds  p  .  Figure  8,  below,  depicts  the 
the  disposition  of  E_  and  E  in  the  (f  ,p)-plane. 


*The  proof  of  the  statements  that  [p]  >  0  implies  \t]   <  0,  and  that  [p]  =0,  if, 
and  only  i:^  [tl   -  0,  will  be  given  at  the  end  of  this  subsection. 


Note  that  dp  =  l-^j     df   +    f  ^j  dS;  we  arrive  at  the  desired  identity  after 


*lote  that  dp  =  I  -^ 

's     v'r 

substituting  for  dS  from  the  basic  relation,  TdS  =  de  +  pdf . 


Figure  8:  Depicting  a  means  of  Joining  a  state  (?l,p^)  behind  a  shock  to  a 
state  (T  »P  )   in  front  by  means  of  an  arc  of  the  constant  internal  energy  curve 
E^  joining  (T^fV^)   to  (^^jp^),  followed  by  the  arc  of  the  (full-lined)  adiabatic 
curve  A^  Joining  (^^jP^)  and  (^-^,p^).  L  is  the  line  through  (?^»P-,)  and  (f  ,p   ), 
The  dashed  curve  A^  is  the  hypothetical  adiabatic  arc  introduced  in  the  proof  of 
case  i)  and  shown  to  be  inadmissible,  while  R.  is  the  curve  of  constant  energy 
through  {T^,p^). 

Now  let 
(3.11)     A^:   p  =  g(r,S^) 


be  the  adiabatic  curve  through  ('^■,,p,).     Since,  from  (3.10),   (3.2)  and  (3.3) 
we  have 


(3.12) 


0     >     (f ) 


e 


<lf'    ' 


we  may  conclude  that  A,  lies  below  E-  as  long  as  ^  >  '^  [see  Figure  8] .  Let 
('^  ,p  )  be  the  point  of  intersection  of  A-  with  E  .  If  ^S  is  the  change  in 


I      t 


entropy  along  the  arc  of  E    which  joins  (  "^ojPq)  ^^  (  ^  'PqX  ^^^^  employing 
the  basic  relation 
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(3.13)     TdS  =  de  +  pd"^ 

we  find,  since  de  vanishes  on  E  .  that 

(3.1U)    /^S  -  (      ^dT. 

h 

This  change  in  entropy  is  the  same  as  the  difference  in  entropy  [s]  of  the 
states  (  ^1  jPt  )  sind  (  2"  ,p  )  since  the  point  (  T'jP' )  is  joined  to  (  Tn  jP-i ) 
by  the  adiabatic  path  A  .  Thus  we  conclude 


(3.15)    [s]  =  f  §dr. 


The  problem  of  proving  [f\  <  0  implies  [sj  >  0  is  therefore  reduced  to 

showing  that  f      exceeds  T  . 
o  o 

That  r'  must  be  greater  than  ?"  is  a  consequence  of  the  Hugoniot 
relation  (3.5)  and  our  assumption  (3*2)  that  p  =  g(T~,S,)  is  convex  downward. 
This  may  be  seen  frcm  the  following  argument:  Suppose  ?'  <  r  .  Then  because 
the  curve  A,  is  convex  downwai^  [see  (3.2)]  it  intersects  the  line  L  through 

(  T  ,p  )  and  (2"  ,p  )  in  at  most  two  points  -  (  ■?"t»P-i)  and  possibly  (  ^jPq)*  I** 

f 
(3.16)    W  »  j   pdr 

be  the  area  between  T   and  ?"'  and  under  A, .  This  area  is  clearly  less  than 

o  o  1 

the  area  between     Ti   and      t"  and  under  the  line  L.     Thus  we  have: 

1      o 


rl 


VPc 


(3.17)   w  =  f  pdr  <  -^  (  r^  -  r^). 


^1 

But,  since  A,  is  an  adiabatic  path,  it  follows  from  (3«13)  that  W  »  ®i~®o  ^^'^ 
hence,  using  (3.17),  that 
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(3.18)  e^-e^  <  ^  (  7-0  -  ^i). 
From  (3.5),  however,  we  find  easily  that 

(3.19)  e^-e^  >  -A^   (  r^  -  r^) 

for  a  compressive  shock.  VJe  have  thus  arrived  at  a  contradiction.  The  assumption 
that  r-  <  7"  is  therefore  false:  t^   exceeds  T^   and  hence  ISJ  >  0  in  a  compressive 

0~0  O  O  *'-' 

shock. 

To  prove  the  direct  part  of  statement  ii)  we  make  use  of  the  fact  that 
the  entropy  may  be  regarded  as  a  function,  S  =  S(e,  2"),  say,  of  e  and  P".     Since 
(e,-e  )  vanishes  when  (  "^t  -  ^)  vanishes  [see   (3.5)]  it  follows  that  the  entropy 
behind  the  shock  is  identical  with  that  in  front.  Thus  [j|  =  0  implies  [sl'O. 
The  proof  of  the  converse  follows  the  proof  of  iii). 

To  prove  statement  iii),  namely  that  \t^  >  0   implies  [s]  <  0,  we 
proceed  as  follows.  Let  us  write 

-^1  -^o'  ^o  =^1 

(3.2C) 

In  terms  of  this  new  notation  we  have  to  prove  that  ['^l  <  0  implies  [sj  >  G. 
The  proof  of  statement  i)  now  applies  to  the  'roofed'  quantities  and  furnishes 
the  desired  result;  namely,  that  an  expansive  shock  implies  a  decrease  in 
entropy  behind  the  shock. 

To  prove  the  converse  of  ii)  we  begin  by  supposing  that  [sJ  »  0  implies 
[r]  ^<  0.  Employing  i)  and  iii),  it  then  follows  that  [s]  ><  0,  contrary  to  our 
original  assumption.  Thus  [s]  =  0  implies  [^1  "  ^* 

l(fe  have  yet  to  show:  l)  [p]  >  0  Implies  [^1  <  Oj  2)   [p]  =  0  if ,  and 
only  if,  L^]  =  0*  Now,  the  statement  that  [p]  >  0  implies  [^]  <  0  is  equivalent 
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to  the  statement  that   [t]  >  0  implies   [p]   <  0.     Since  we  have  found  that 
SL   <  S     for  an  expansive  shock  it  follows  directly  from  (3.2)  that 
p^=  gC^r^,  S^)  <  gC-t^,  S^)    =  p^  whenever    f^  >    t^j  this  proves  1). 
Let   [tj    =  Oj  then  we  know  that   [3]  =  0.     Consequently, 

Let   [p]   =  0  and  suppose  that   [f]  ><  0.     It  follows  that  [s]  ><  0. 
But  this  result,  using  (3.2)  again,   implies    [p]  ^<  0,   contrary  to  our  original 
assinnption.     Thus   [p]    =  0  implies   [p]    =0,  which,  together  with  the  previous 
result,  proves  2). 

3,2  The  relation  between  d[s]   and  dm  in  a  compressive  shock 
The  object  of  this  subsection  is  to  prove   : 


[^Hl 


2     1 


(3.21)  T^d[s]   =     <    [tY  *  —2 


n 


"    mdm,  0  <  6  <  e-| . 


This  relation  is  valid  for  all  compressive  shocksj  the  second  term  in  the  braces 
vanishes  when  H  vanishes  or  when  the  shock  is  a  pure  gas  shock. 


We  begin  our  proof  by  observing  that  the  result  of  eliminating  [u  ] 


Again  we  are  asstming  that  the  transverse  magnetic  field  H,  in  front  of  and 
behind  the  shock  is  parallel  to  the  y-axis.  This  assiimption  is  justified  in 
Section  li. 
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from       (1.3  -h^  „)by  means  of  (1,3-A-)     is 


(3.22)  n\t\  ♦    [p  +  nH W2]     -    0. 

•7 


If  we  eliminate  [u^_]  for  (1.3  -  A^  .)   by  means  of  (1,3  -  Ao  t  ^ 
we  find,  in  addition,  that 

(3.23)     m^LtH^l  -  HH^  [Hy]   =  0. 

Taking  the  differential  of  (3.22),  (3.23)  and  of  the  hydromagnetic 
Hugoniot  relation  (3.5),  we  obtain  the  following  differential  relations  s 

(3.2U)      0  =  dm2[t]^  -  [p  +  tiH^M  d[rj  +  M  d[p  +  piH^  /2I  , 

(3.25)  0  =  dm^r/Hyl^y^^     ^  ^|-j^^-[   ^j-^y   _  ^[-^y  ^j-j^^j    ^ 

(3.26)  0  =  d[e]  +  ?  d[t]  -^   [r]   d[p]/2  +  h[h j2  d[tr]A  +  iid(  [h J^)  [r]  A. 
Employing  the  relations, 

(3.27)  T^  dS^     -     Ttl  d[s]    =  de^  +  p^d  ^^  =  d[e]   +  P]^d[r], 

(3.28)  [H^      -     [Hy]2  .  2Hy^^    [h^]. 


In  the  following,  A.       and  A.    .    ,   j  =1,2,  are  the  normal  and  transverse 
component  equations,   respectively,  of  A.. 

When  H     =  0  or  when  the  shock  is  a  pure  gas  shock,  each  term  in  (3«23)  is  zero. 
In  these  cases,   (3.25)  reduces  to  0  =0, 
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equations  (3.26)  and  (3.2U)  become 

(3.29)         2T^d[s]  =  dW  ][p]  -  h[h  ]V2[    -  [r]  d[p]  -  ^(M/2)d([HJ2)  ^ 


(3.30)       Ity  dm^  -  d[r] 


M   -  4Hy]V2[-    [tjdfp]  -  j.W/2  d(^yl2)  -  ^iMHy^^d[Hyl 


while  (3.25)  becomes 

(3.31)         dm2    [fHylVnH^-   =     Hy^JfJ  d[ig   -  Hy^^   [h^]   d^. 

If  we  subtract  (3.31)  and  (3.p)  from  (3*29),  we  obtain  the  desired  result, 
namely  (3,21), 

U,  Magnetic  shocks  -  preliminary  remarks 

Magnetic  shocks,  it  will  be  recalled,  are  compressive  shocks  in  which 
|_H.  ]  does  not  vanish.  In  this  section  we  shall  rewrite  the  shock  relations 
(1.3)  in  a  form  most  appropriate  for  dealing  with  magnetic  shocks. 

We  begin  by  showing  that  there  is  no  essential  loss  of  generality  in 
assuming  that  H.   and  u .  both  in  front  of  and  behind  the  shock  are  parallel  to 
the  y-axis.  Toward  this  end  we  eliminate  [u+_]  from  (1.3  --A-,)  by  means  of  the 
transverse  part  of  (1,3-Ap)  -  that  is,  by  means  of 

(U.l)     m[u^J  .  ^^\L«t^"l 

to  get,   assvaning  H     /  0, 

(U.2)  (m^Ti  -  ^xl^)  H,^,i  »  (m^r^  -  ^^)  H,^^^,        T^  <   T  . 
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We  consider  two  cases:  case  a),  H.  „  •»  0:  case  b)  H.    /  0,  In  case  a) 

T;r,o    '         tr,o 

the  factor  m  T  -  llh  must  vanish:  otherwise  H.   -  would  have  to  vanish  for 
1    n  "tr,x 

aU.  values  of  the  shock  strength  parameter  and  we  would  not  have  a  magnetic 

shock,  /.ssuming,  then,  that  m  f  -  jjJT  vanishes,  it  is  easy  to  verify  that 

the  shock  relations  detentiine  only  the  magnitude  of  H.   ^ ,  leaving  the  direction 

tr,  j- 

arbitrary.  We  are  therefore  at  liberty  to  choose  H.  ^  parallel  to  the  y-axis« 

xr,  J. 

In  case  (b),  it  follows  directly  from  (Ij..2)'^  that  H.   ,  is  collinear  with  H^ 

tr,x  tr,o 

If  we  take  H.    parallel  to  the  y-axis,  H. ^  ,  will  also  have  this  property. 
Assiiming  that  this  choice  has  been  made  it  follows  frcsn  (U.l)  that  u  ,  -  u 

Z,l     Z,0 

vanishes.  By  choosing  a  suitable  (moving)  coordinate  system,  it  is  clear  that 


u  ^  and  u    may  be  assumed  to  vanish. 


We  have  thus  far  supposed  that  H  ^0.  If  H  vanishes  we  find, 
employing  (1.3-.'  )  and  (1,3  -  ■^o  t  ^*  *^* 

(^•3)    '^Ar,l  =  '^o  «tr,o 


*Since     ^^  <    '?^  the  factors  m^  f^-  nH^  and  m^f^-  nH^  in  (h.2)   cannot  vanish 

2  2 

simultaneously.     If  m    T^-  [jE^ji  0,   then  the  right-hand  side  and  hence  the 

left-hand  side  of  (U.2)   does  not  vanish.     Thus  m^"^-  uH^  /  0  and  H^     ,    = 

1     ^  n  '^  tr,l 

\t  o^'^^^o"  t^f^n^Aro^*!-  1^1^).     If»  on  the  other  hand,  m^?^-  nh^^  =  0  then 
ra  '^"  ^^  ^  ^  ^^^  hence  H       ^  must  vanish;  here,  too,  H       ,  may  be  regarded 
as  collinear  with  H^^  ^.     Actually,   as  we  shall  see  later,  m  T-  ^iir  vanishes 
only  for  isolated  values  of  the  shock  strength  parameter. 
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(U.5) 


Clearly,  there  is  again  no  real  loss  of  generality  in  assuming  that  H.      . 

and  u.      jj   j  =  0,1  are  parallel  to  the  y-axis. 
tr,  0 

If  we  now  set 

(U.U)  h     =     (Hy^^  -  Hy^^)/H^, 

and  make  use  of  the  above  result,  the  shock  relations  (1.3  ^-I'^o)  reduce,  using 
the  notation  introduced  in  (l.U)  to  the  following  set  of  relations: 

^2  n   m^^"!^  -  t^(X+  h  sin©^), 

Bg     m^^^(h  -  7i   sin  9^)   =  |iH^  cos^^  h, 

B3    [u  J  -  -  m  ^3^  ^  . 

To  complete  this  set  of  equations  we  must  add  (1. 3-Ai  )  or  an  equivalent  relation. 
Now  it  is  knowrr  -I  that  (1.3-^-, -A.)  and  A.  imply  the  Hugoniot  relation  (3.5) 
when  m  /  0.  By  retracing  the  steps  of  this  derivation  it  is  easy  to  verify 
that  (1,3-A^-A-)  and  the  Hugoniot  relation  imply  (1.3-Ai  ).  The  Hugoniot 
relation  (3,^)   may  therefore  be  used  in  place  of  (1,3-Ai  ),  If  we  emplcythe 
(polytopic)  ideal  gas  relations  (1.3'),  the  Hugoniot  relation  becomes,  in  the 
present  notation: 

2  '?  s^         .2 

(U.5)    B.     X  -  2_  +  S . 

^  2-(y-1)^    2-(y-1)  1^ 

Note  that  (U.5  -  Bp  -  Bi  )  are  invariant  under  all  changes  of  sign  of  h 

and  of  ©  (©  /  0)  which  leave  the  relations  h  sin  B^  >  0   and  h  sin  9^  <  0 
00  0  o 

unchanged,  while  the  expression  for  [u  ]  merely  changes  sign  with  h.  In 

analyzing  M„  and  M  shocks  it  is  clear  that  we  may,  without  incurring  any 
X       s 
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essential  loss  of  generality,  restrict  the  ranges  of  h  and  €•  as  follows 

(U.6)     M^     h  >  0,   0  <»^  <  90*^, 
(U.7)     M      h<0,   0<9   <  S0°. 

In  the  next  two  sections  we  shall  regard  (U.6)  and  (U.7)  as  the 

definitions  of  M-  and  M     shocks  respectively. 
I  s 


5.     Fast  magnetic  shocks  (h»  >  0,  0  <  ©^  <  90°)j  M-  shocks 

5»1  Fast  magnetic  shocks  of  types  1  and  2;  Mj.       and  M^       shocks 
Our  analysis  of  M.  shocks  is  based  on  the  following  relations: 

(5.1)  Wf_-    1]^_^ If bj^./p^^, 


(5.U) 


(5.5) 


Wf         ^n.l     ^f 


h. 


(5-^)      ff^  -  i^  -h-  -  linr  t  -  (7f^)^^  «or''. 


ll/2 


b     ,  V     ,      cos  &  o 

n,l  n,l  o 


^f     k2  /yf    \2  ^2 

(5.3)       (-^1      -    ^U-°^      -  ^ 


^^n,i;  ^f  V"'°/  (P^n,l>'         [l-(^fAf )3in  B,]  ' 

^f /hf  -  3ln  e^       ^     X^  ^  tjh^sinyd-e^)  1/^0^]. 
l-(»[^/h^)sin  e^          ^f  2  sin  ©^-(y-Dh^ 

X^/h^  +  sin  e^  )(^  [  -  J  h^sin  ©o-(l-s^)  1  /^HJ^  ] 

l+(Zj./hj)sin  e^  f  2  s^sin  6Q-(Y-l)hj 
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(5.6) 


d(X^/h^)         (Y-l)(X^/h^)^  +  Y  sin  9^(X^/h^)  +  1 


dh. 


t2fT;(Kp 


(5.7) 


2 

cos  e 


d(;rfAf)  (l-(T^/h^)sin  9^)' 


where 

(5.8)  a)     R     (h^)  -  fe  hj^sin  »^  +  (l-3^j)P  +  (2s^sin  e^-(Y-l)h^)(h^+2  sin  9^)  , 

b)     ^h^)  =     %  V^n»o  "  ^^'^0^1^  *  ^^  sin&^-(Y-l)hj)(h^+2  s^sin  9^). 


(5.9) 


It  is  easy   to  verily  that 


R_(h^)  =  I^(h^) 


R(h^) 


=  h' 


2       2 
Y  sin  © 


-  (y-1) 


+  hjsin  9^(2-y)(1+3^) 


(1-s^)     +  Us^sin^^. 


Equations  (5.l)-(5.3)  follow  almost  immediately  from  (U.5-B^  „,   B- 

and  B-  )  and  the  fact  that  m  =  p-V  -,  =  p,  7„  .  The  first  equation  of  (5«U) 
2,y  1  nji    o  n,o 

is  obtained  by  eliminating  ra  t^   from  •(U,5-Bo  _)  by  means  of  (U.5-B«  ),  The 
first  equation  of  (5.5)  is  obtained  from  this  result  simply  by  solving  for 
Jf-Ai^  in  terms  of  X„/h„.  The  last  equations  in  {S,\x)   and  (5.5)  are  obtained  as 
follows t  When  ^./^^   is  eliminated  from  the  first  equation  in  (5.U)  by  means 
of  the  Hugoniot  relation  (U»5-Bi  ),  it  is  found  that 

.2  '' 


(5.10)   l^ 


2  sin  »^  -  (Y-l)hj 


2(l-s^)-YhjSin  9^[  -  (h^+2s^sin  9q)-0. 


Similarly,  employing  the  Hugoniot  relation  to  eliminate  X«/h-  from  the  first 
equation  of  (5.5)  it  is  found  that 
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(;.n,  (|A'! 


2  s^sin  »^-(Y-l)hf 


>•(!;: 


Yh^in  e^+ 2(1-3^) 


-(h^+2  sin  0jj)=O. 


Vw 


The  desired  results  -  the  last  equations  in  (5.U)  and  (5.5)  -  are  obtained  by- 
solving  these  equations  for  X_/h»  and  ^^/h^.     To  derive  (5«6)  we  differentiate 
(5.10)  with  respect  to  h.,   solve  for  d(X-/h-)/dh-  and  employ  the  last  equation 
of  (5.U)  to  remove  X-/h  frOTi  the  denominator  of  the  resulting  expression. 
Equation  (5.7)  follows  immediately  from  the  first  equation  of  (5.U)  on 
differentiating  X^AfWith  respect  to  t'^Ai^.. 

The  expressions  for  >7f  Ax.  and  X-A„  which  result  from  choosing  the  plus 
branch  before  the  radical  in  (,$.h)  and  (,^,^)  will  be  referred  to  as  plus 
branches     and  will  be  denoted  by  Xj./h„  and  <\^/'h',  the  remaining  branches  will 
be  called  minus  branches,  and  denoted  by  X^/h^  and  »f^Af.     Employing  (5.6)  and 
(5.7)  we  obtain  the  following  useful  result:   Th£  h^-derivatives  of  X/h^  and 
>f „/h„  are  positive  or  negative  according  as  we  are  dealing  with  the  plus  or 
minus  branches}   it  is  assumed  here  that  h_  is  limited  to  those  values  for  which 
/U  is  real, 

Equations  (5.1)-(5.3)  and  the  last  equations  in  (5.U)  and  (5.5) 
determine  a  one-parameter  family  of  states  behind  the  shocks,  with  h„  playing 
the  role  of  the  shock  strength  parameter. 

We  turn  now  to  the  problem  of  ascertaining  the  admissible  range  of  h_, 
that  is,   the  totality  of  values  of  h-  for  which  the  above  one-parameter  family 
is  an  admissible  (see  Section  1.2)  solution  of  the  shock  relations.     In  light  of  our 
definition  of  admissibility  and  the  results  of  Section  3,  it  is  clearly 
sufficient  to  determine  the  range  (or  ranges)  of  h_  in  which  a)  the  state 
behind  the  shock  is  real  valued j  b)  the  shock  is  compressive,  i.e.,  »^«  >  Oj 

itate  behind  the  shock  depends  continuously  on  h„  and  -t 

Conditions  a)  and  b)   are  equivalent  to  the  conditions 


c)  the  state  behind  the  shock  depends  continuously  on  h.  and  the  state  in  front. 


*It  is  easy  to  verify,  using  the  first  equation  of  (5.U),  that  >^/h^  corresponds  to 
X»/h-  and  "^/^r.  corresponds  to  Xl/h„. 
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a')   R(h^)  >  0 

(5.12) 

b«)   sin  9^  <  *LAi^  <  l/sin9^. 

Condition  b')  is,  in  turn,  equivalent  to  the  relation 

(5.13)    b")    Xj/h^  >  0. 

To  prove  these  equivalences  we  observe  first  that  if  X-/h„  and 
)^„A„  are  real,  then  a')  is  valid.  Second,  we  note  that  X-/h_  (=  Y^^SQAf.  +  hV2) 
is  positive J  for,  h^  >  0,  and  "t _  >  0  imply  Y _  >  0 .  Employing  this  result  and 
the  first  equation  of  (5.1i)  we  find  that  X„/Ai„  >  0  if,  and  only  if,  the  relation 
(5.12-b')  holds.  Condition  a')  and  b')  are  therefore  necessary.  That  a')  and 
b')  are  sufficient  may  be  concluded  from  a  siirilar  argument.  Below,  we  shall 
employ  a'),b'),  or  b") and  statement  c)  above  to  deteimine  the  admissible  range  of  h  , 
We  begin  by  rewriting  the  second  equation  of  (5.U)  as  follows: 


(5.11;)     r^ 


^f    ^1*'/^ 


^f       °X 


where 

(5.15^     4  "  7  ^f  ^^  ®o  "  ^^"^o^  ' 

(5.16)  °x  *=  2  sin  »^  -  (Y-l)h^  , 

(5.17)  ^  =  ^X  "^  (2  sin&^  -(Y-l)h^)(h^  +  2  SpSine^). 

Now  let  us  restrict  h-  to  the  range 

(5.18)  0  <  h^  <  h^  =  2sln9y(Y-l). 

In  this  range  the  product  of  the  last  two  tenns  in  the  right-hand  side  of  (5.17) 
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is  positive^  hence  /rI     exceeds  |byj.  In  addition,  Cj.   is  positive.  Clearly 
then,  only  the  plus  branch  leads  to  positive  values  fori  /h  and  conversely. 
Since  '?f./h^  corresponds  to  X./h  [see  last  footnote},  it  follows  that  only  the 
plus  branches  are  admissible  in  the  range  0  <  h„  <  h-. 

We  turn  now  to  the  problem  of  determining  under  what  circtmistances  it 
is  possible  to  continue  the  plus  branch  into  the  range 


(5.19) 


^f  ^  ^f 


2  sin  9  Ay-D  . 


To  this  end  let  us  evaluate  the  function  i7„/h„  at  h«  =  h„5  it  is  easily  verified 
if  we  use  the  second  equation  of  (5.5),  (5.9),  and  (5.17)>  that 


-  + 

(5.20)    r^ 
f 


¥^f 


l/sin  ©, 


Y  sin  e. 


if 


A 

h. 


>o, 


(y-1)(1-s J  '    ^^    ^ 


h^=h^ 


<  0. 


The  curves  passing  through  (l/sin  ©^)  at  h^=  h^  will  be  referred  to  as  type  1 
curves,  and  those  passing  through  y  sin  9^/(y-1)(1-3q)  at  h^=  h.  as  type  2  curves. 
These  fonn  two  mutually  exclusive  families  characterized  by  the  relations 


(5.21)     type  1:    s^  >  1  -  y  sin^^/(Y-l), 
type  2:    s   <  1  -  Y  sin  9^/(y-1)> 


these  are  simply  the  conditions  b^ 


>  0  and  b- 

V^f " 


<  0  written  out  • 


hf=h. 


Since  X_/h  becomes  infinite  when  *tf/^f  approaches  l/sin  &  Lsee  (5.U)J, 

A 

it  is  clear  that  type  1  curves  cannot  be  continued  into  the  region  h„  >  h-.. 
For  type  2  curves,  y sin  ©  /(y-1)(1-s  )  [aee  (5.20)]  is  clearly  positive; 
(5«12-b')  is  therefore  automatically  satisfied.  Consequently,  using  (5.U) 


-  66  - 


again,  X^A^.  is  finite  and  positive.  In  addition,  it  follows  from  (5.6)  and 
(5.7)  that  d(4'-/h-)/dh„  is  a  finite  positive  value.  Clearly,  then  type  2  curves 
may  be  continued  into  the  region  h.  >  n-. 

Let  us  determine  an  upper  bound  for  this  range  by  finding  the  inaxiiiiuia 
value  for  which  R^  >  0,  i.e.,  for  which 

(5.22)     b^  -  ((Y-l)hj  -  2  sin  9o)(h^  +  2sQ3in  9^)  >  0. 

To  this  end,  we  observe  that  in  curves  of  tjrpe  2,  b  is  zero  at  a  value  of 

h-,  h^,  say,  which  exceeds  h-j  for  b  vanishes  at  the  value  h-  =  2(1-3jj)/y  sin  9^ 

which  clearly  exceeds  h„=  2  sin  9  /(y-1)  when  s^  <  1  -  y  sin  9^/(y-1) 

A 

[cf.  (5«21),  type  2].  In  addition,  we  note  that  Ry  is  positive  at  h-  =  h_ 
and  decreases  monotonically  to  a  negative  value  at  h«  =  h-j  for  b„  decreases 
while  the  two  factors  in  the  last  term  of  (5.22)  decrease  as  h«  increases  from 

A         I 

h_  to  h-.  Thus  Ry  has  a  zero  between  tliese  limits.  Now,  as  a  mere  inspection 
of  the  coefficients  of  (5.9)  will  reveal,  IL.  has  a  positive  root  if,  and  only  if, 

sin^9^  <  U(Y-l)/y'^. 

Setting  the  right-hand  side  of  (5.9)  equal  to  zero  we  find  the  following 
expression  for  this  positive  root: 


A    sin  9  (2-r)(l+3„)  +  2  cos  9^  /(y-D ( 1-s^ ) ^+  s  Y^sin^ 

(5.23)  \ 2 2 o^ 2 2 £ 

^  2(y-1)  -  (y  sin\/2) 

Thus  ^R     is  real  in  the  range  ti^  <  h-  <  h,.  . 


This  relationship  is  satisfied  by  a  type  2  curve  for  values  of  y  <  U. 
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Two  simple  facts  are  decisive  in  determining  whether,  in  curves  of 

^        <%  _    .   , 

type  2,  the  entire  range  h-  <  h_  <  h„  is  admissible.  The  first  is:  /K-  <  |by|, 

h  <  h „  <  h-.  This  insult  follows  immediately  from  (5.22)  because  one  factor 

in  the  product  term  is  positive  while  the  other  is  negative.  The  second  Is: 

b_  <  0,  h „  <  h _  <  h~.  This  is  immediate  for  curves  of  type  2j  in  fact,  b-  is  negative 

in  the  larger  range  ^f  5  h _  <  h-,  where,  it  will  be  recalled,  h„  is  the  zero  of  by. 

Now  let  us  consider  the  expression  for  X_/h-  in  (5.1h).  In  the  range 
h _  <  h -  <  h-.,  Cj  is  negative  and,  as  we  have  just  seen,  jKy.  <   |by|.  Consequently, 
using  (5,lli),  we  find  that  X  /h-  >  0  in  the  given  range  if,  and  only  if,  by  <  0. 

It  follows  that  X-/h-  is  admissible  in  the  entire  range  h-  <  h«  <  h..  Moreover,  both 

plus  and  minus  branches  of  X„/h-  and  hence  of  'T^/h^  are  admissible.  X„/h-  and 

—  <^       'ft 

vf„/h„  are  thus  double-valued  functions  in  the  range  h„  <  h .  <  h,.,  their  minus 

branches  joining  the  corresponding  plus  branches  at  endpoint  of  the  range, 
namely,  at  h.  «»  h„. 

5,2  Description  of  the  curves  X_/Ai.  versus  h_  and  i^/^f  versus  h_ 

5.2.1  M^-'-^  curves  (s^  >  1  -  y  sin^e^Ay-l) ) 

The  solid  line  curves  of  Figures  9a  and  9b  depict  curves  of 

t3rpe  1.  .'^s  noted  above, type  1  curves  exist  only  in  the  range, 

^     2  sin  e 
(5.2U)     0<h^<h^  =  -__£  , 

and  within  this  range  only  the  plus  branches  are  admissible,  X-A^-  assumes 
its  minimum  value  at  hf=  0,  namely, 

/(l-3^)2  +  Us^3in\  -  (1-s^) 


4 

(5.25)    T^i 
*^f 


k  =n  2  sin  ©„ 

h  =0  o 


-  68  - 


\/\  vs  hg 


1(l-SQ)^+4s„sin2e_-(l-s„) 


0  0  ° 


2  s„  sin  6^ 
0  o 


h^=  2  sin  B^/iy—l) 

^  ZcosSq    V(y_|)(|_s^)2+   y2s^sineQ+    (2- y)(l +So)sin  5^ 


^     = 


2(y-l)  —  y^sin^%/2 


X  /h     vs  h 


So>l-rsin  9o/{r-\) 


^  Figure  9(b) 


2  s  i  n  fl. 


s^<\-ys\n^8^/(y-\) 


•[l+Sp(y-l)]sine^ 
(l-So)(/-")-rsin%o 


I 


A 

h 
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FJgraires  9(a),  9(b).  Sketch  of  the  variation  of  ^^Ai^  with  h^  (a),  and  of 
X„/h-  with  h_,  (b)  in  M^  shocks.  The  solid  line  curves  are  of  type  Ij  the 
broken  line  curves,  type  2.  The  (±)  superscripts  indicate  the  plus  and 
minus  branches.  X-/h„  and  "^^A^r.  "^ary  in  the  same  sense  as  h^  on  all  plus 
branches  and  in  the  opposite  sense  to  h-  on  all  minus  branches. 
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and  increases  monotonically  to  -»-oo  [see  (5.6),  plus  branch],  approaching 
the  line  h^-  h„  asymptotically  as  h-  increases  to  h».  The  ratio  IfA- 
varies  in  the  same  sense  as  X_Ai,  [see  (5.7)] J  it  increases  monotonically  from 
the  value 


(5.26)  ^ 


yH^sf  *  Us^sin\     -  (1-s^) 
h^O     '  2  s^sin  d^ 


to  its  maximum  value  (l/sin  9  )  at  h-  =  h„  as  h-  increases  from  zero  to  h». 
Note  when  h-  =  h«  that  »2-  ■  2/(y-1)  which  is  the  same    maxinnira    value  achieved 
in  pure  gas  shocks. 

5.2.2  M^^^  curves     (s^  <  1  -  y  sin^e^Ay-l)  ^ 
Curves  of  tjrpe  2  exist  in  the  range 

(5.27)  0     <     h^    <    h^ 

where  h-  is  given  in  (5.23).  In  the  range 

(5.27')    0  <  h^  <  h^  =  2  sin&y(Y-l) 


^ 


only  the  plus  branches  are  admissible;  in  the  range  h„  <  h»  <  h«  both  plus 
and  minus  branches  are  admissible.  Since  c-  <  0  in  the  range  h_  <  h.  <  h», 

the  values  (X-/h»,  ^-~/h„)  exceed  the  corresponding  values  on  the  plus  branch 

^  +  _+ 

except  at  h_  »  h»  where  the  branches  join.  Since  d(X^/hJ/dh-  and  d(^A'^)/dh- 

are  all  infinite  at  h-  =  h-,  it  is  clear  that  this  joining  is  a  smooth  one. 

_,      .  <% 

X_/h-  and  l^'/h-  assume  minimum  values  at  h_  «=  h»  and  inci^ase  monotonically 

[see  (5.6)  and  (5*7)1  to  -t-oo  and  (l/sin  ©  ),  respectively,  as  h-  decreases  from  h_to 

h-.  On  the  plus  branches «  ^^/h-  and  X-/h-  assume  minimum  values  at  h-  =  0 

[see  (5.25)  and  (5.26)]  and  increase  montonically  [see  i$t6)   and  (5.7)^  to 

maximum  values  at  h-  »  h_.  Note  that  the  zeros  of  the  denominators  in  the 
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expressions  for  >(-  A  and  I-A  occur  at  the  zeros  of  the  corresponding 
numerators.  These  zeros  therefore  do  not  produce  infanities:  in  fact, 
employing  L'hospital's  rule  or  solving  directly  from  (5.10)  and  (5.1l)»  we 
find  that 


(5.28) 


b) 


^ 


sin  &„[8„/ff-D  +  d 


o*-  o 


V  ^  ^^^°  ®o    ^-^o  *  ^/y-1  s^sin^e^ 


hf  =  Y=I  ^^"  ^0 


1-3  -(Y/r-i)'3in^^ 


The  dashed  curves  of  Figures  9(a)  and  9(b)  depict  typical  curves  of  type  2. 


5,3  Some  general  features  of  M-  shocks 

Our  chief  object  here  is  i)  to  prove  that  Jf-,  Y„,  v  ,/h     -,  V  ^/b„  ^, 
"^  f   1   n,l  n,l   n,o  n,o' 

m»  and  [s]-  vary  in  the  same  sense,  and  ii)  to  determine  how  these  quantities 
vary  with  h».  In  view  of  (5.3)  and  (3.21),  i)  will  be  proved  if  we  can  show 
that  >f „  varies  in  the  same  sense  as  ^»/h»  and  Y-  varies  in  the  same  sense  as  Tf 

Consider  first  the  variation  of  yf„  with  'L/bf •  ^s  we  showed  in  the 
preceding  subsection,  "^^h^  varies  directly  with  h-j  this  statement  is  independent 
of  type.  Thus  >it  =  (^/h-.)h-  varies  directly  with  ?t,/h-.  It  remains  for  us  to 
show  that  the  same  is  true  for  the  minus  branch.  Our  procedure  will  be  to 
prove  that  d'^l/dh  has  the  same  sign  as  d(if"Ai£.)/dh„,  which,  in  virtue  of 
(5.6)  and  (5.7)  is  nei^ative.  Now,  evidently  we  have 


1  '^'^f 


""dhl 
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since  d(»f"/h„)/dh-  is  a  large  negative  number  in  the  vicinity  of  h-  «»  h-, 

it  follows  that  dlj»/dh-  is  also  a  large  negative  number  in  this  vicinity. 

To  obtain  the  desired  result,  it  is  sufficient  to  prove  that  7fAf  exceeds  (sin9  )" 

when  d^r/dh-  vanishes;  for  then  dTzT/dh-  would  be  negative  through  the  admissible 

range.  To  this  end  we  observe  that 

when  d*^»/dh-  vanishes.  Differentiating  (5.11)  with  respect  to  h-  and  eliminating 
d{n^/h  )/dh^  by  means  of  (5.28),  we  find 


(5.31)     (f)    [Ls„  Sin  V  (r-Dh,]  =  -  2(1-.)  (|)  -  "f 


The  addition  of  the  left-hand  side  of  (5.11)  to  the  right-hand  side  of  (5.31) 
leads  directly  to 

The  requirement  that  the  right-hand  side  of  (5.32)  be  positive  results  in  the 
desired  relation,  namely. 


S  =■  r-nrw:  i^^-rr-)  >  (s-V'- 


sin  e. 
..^    Y  sin  9^  ^"  ■   hj 


Our  aim  now  is  to  prove  that  Y-  varies  directly  with  vf  _.  As  an 
inspection  of  the  Hugoniot  relation  (U,5  -  Bi  )  will  reveal,  T.  increases  with 
increasing  4 _  as  long  as  we  restrict  ourselves  to  the  plus  branches  of  T„  and 
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_+ 


f^-j  this  is  so  because  >^«  varies  directly  with  h-.  On  the  minus  branches  of 

(2) 
M«  '  shocks  the  behavior  of  the  second  term  is  no  longer  decidable  by  a  mere 

inspection^  for  here  >t-.  and  h«  vary  in  opposite  senses.  Nevertheless,  Y_  varies 

directly  with  ?^  on  this  branch  also,  as  we  shall  now  show.  Firsts  we  observe 


that  the  IZ.  derivative  of  XZ     is,  using  the  first  equation  of  (5.U), 

|l-  ^  sin  »  I   l-sin  »Q/(d"^-/dh^)|  +  fe"  -  h^sin  »Jsin  &gd(?-/h^)/d"^- 

■       B      •  _____ 


(5.33)      -^ 


f  1  -   (^-A^)sin  9   I    . 


The  right-hand  side  of  this  equation  is  positive,   for  fein  ^q)"    >  ^f/'^f  ^  ^^°  ®'o 
[cf.   (5.12-b')]  and,   as  we  have  seen  above,  d(*f*/dh_)  and  d(»J'/h„)/dh-  aire  both 
negative.     The  relationship  between  the  derivatives  d?"/d"^'  and  dX»/d»?*  may 
readily  be  obtained  by  differentiating  (l,Li-e).     The  result  is  found  to  be 


(5.3U)   -^  =  ^ 


An}      ^^ 


dX^     h^ 


dJ?"     drf^/dh^ 


Since  dXl/dW"  is  positive  and  d'^/dh.  is  negative,  df'/dfX  exceeds  zero,  which 

is  the  desired  result. 

It  is  now  a  relatively  simple  matter  to  prove  ii),  that  is,  to 

determine  how  ff-,  T»,  etc.,  vary  with  h-.  To  this  end  we  show  that 

a)  if-,  varies  with  h„  in  the  same  way  that  ^-/h-  varies  with  h-  [see  previous 

subsection  and  Figure  9^)1  and  b)  that  Y.,  V^  ^^  t,  V^  ^^„  ^,  m„  and  [oj- 

-^        ~*— ""  I   ^y*^  **5-^   n,o  n,o   X  ^■^■*    i 

vary  with  h„  in  the  same  way  that  X»Af.  varies  with  h„  [see  Figure  ^tj;^ .  Note 
that  statement  a)  involves  variables  of  finite  range,  statement  b)  of  infinite 
range;  the  variation  with  h-  of  both  types  of  quantities  is  the  same  except  for 
this  fact. 

To  prove  a)  we  merely  recall  that  5[-  varies  in  the  same  sense  as 
ilj/h».  To  prove  b)  we  note  that  d(X^/h^)/d(  ^^/h^)   exceeds  zero  [cf .  ($.^jj 
so  that  Y„  ,  V  n  A>_  n  etc.  vary  directly  with  X^/h». 
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Finally,  we  would  like  to  show  that  ^  and  Y„  are  independent  of  & 
when  h«  =  2(l-s  Yf ,   s^  <  1.  To  this  end  we  note  ^T^/hj  =  1  when  X-./h-  =  1 
[see  (5.5)3.  Substituting  for  h„  or  VJ.  and  X-  in  (li.5-Bi  )  and  solving  for  h- 
we  find  that 


(5.35)     h^  -  2(1-s^)/y  , 


«o  <  ^' 


(5.36)    ff, 


2(l-s^)/f, 


So<l, 


and,  using  (l.li-e),  that 


(5.37) 


(y  +  3^  -  1),     Sq  <  1. 


5»li  Weak  fast  shocks 


— + 


If  the  expression  for  ^  in  {.S.^)   is  expanded  in  the  neighborhood  of 

f 
li„=  0,  it  is  found,  after  carrying  out  suitable  algebraic  manipulations,  that 


(^•38)      ^f  =  Ss^f*^ 


o      "o 


(Y-l)(l-r^)   (y-1)(1+Sq)  -  Ys^r^ 


sin© 


1  +  s^  -  2s^r» 
o    of 


h  „  +  . .  • 


where 


(5.39) 


1   1  *  ^0  -^ 


/(1+s^)^-  Us^cos^e^ 


2  cos 


^9 


"f.O 


n,0 


>  1, 


[cf.  (1.5)]. 

A  sufficient  condition  for  the  convergence  of  the  above  expansion  is 
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h  h 

(5.kO)  <   -^  ,  ^ 


iHUS 


"      /o^V^T^V^y 


<<  1,  &  >  0. 
'    o 


If  we  invert  the  power  series  above  to  obtain  h-  as  a  function  of  ^»  and 

substitute  the  resulting  expression  for  h„  into  equations  B,   -  B_  of  (li.5), 

we  find,  after  expanding  with  respect  to  >^_  and  retaining  only  the  leading  terms, 

that 

sin9. 


(5.1'1) 


"=1 

''f' 

T. 

r^       I 

••  i 

\7 

[".If 

= 

r^c^yan 
1-r, 

^0 

•2f 

+  . 

'•  } 

"2 

"f^l 

= 

f,0 

+  . 

•  • » 

S 

Hf 

a 

-°f  ,0  '^'f 

+ 

•  •  • 

» 

0 

<%< 

50°. 

In  addition,  C^  and  (2,70)  yield  the  relation 

which  is  the  same,  up  to  terms  of  the  second  order,  as  the  corresponding 
expansion  in  a  gas  shock.  The  fourth  term  in  (5.1'2)  represents  the  lowest 
order  contribution  axlsing  from  the  presence  of  the  magnetic  field. 

It  is  well  known  that  the  increase  in  entropy  behind  a  pure  gas  shock 
is  of  the  third  order  in  the  shock  strength  parameter.  A  similar  statement 
may  be  made  for  a  fast  magnetic  shock.  Now  the  increase  in  entropy  of  an 
ideal  polytropic  gas  which  passes  from  the  state  (  ?■  ,p  )  to  the  state  (  '^■i  >Pn ) 
is  given  by 
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(5.ii3)  [s]. 


Cy    log 


Pi 


Pl\-Y 


c^  \  log(Y^  ♦  1)  -  Y  log(lf^  +  1) 


Substituting  (5.U2)   into  the  first  term  in  (5.U3)  and  expanding  the  resulting 
expression  in  the  neighborhood  of  7„  =  0,  we  find 


y(y-i) 


(5.UW  [s]^  =  c^  ti^ 


sin^9„ 
Y-*-l  o 


~^  ^  (l-rj^s 


f'  "o 


ifj.    0  0      ,         9^>0. 


The  second  teim  in  the  brackets  represents  the  contribution  of  the 
magnetic  fieldj  when  this  term  is  absent  (5. hi;)   reduces  to  the  corresponding 
expression  for  a  weak  pure  gas  shock. 

5,5    Relation  of  the  normal  flow  velocity  (relativB  to  the  shock)   in  front  oi  and 
behind  fast  shocks  to  the  disturbance  speeds  in  these  regions 


Employing  (5.3)  and  the  second  equation  of  (5.5),  we  find 

1+s  +  /(1+s   )     -  Us  cos  © 
0    /  o  o  c 

jn       V  =0 

f  '  f 


(5.1x5)  lijn      vf  ,      -     liin       vf  ^  =  b^  ^  ■ 

h,->0     "»^         h.->0     "'°         "»° 


[1/2 


2  co^» 


Comparing  the  last  expression  on  the  right  with  (1.5)  we  find  that  it  is 
identical  with  the  fast  disturbance  speed  c„       in  front  of  the  shock -i.e., 

1,0 

(5.1-6)     lim  V^  ,  =  lira   V^   =  c.  , 
h^->0  "'^   h^->0  "*°    ^'° 

In  subsection  (5.3)  we  saw  that  V  /b         increased  with  increasing  w»:  clearlv 

n,o  n,o  ^  'f 


(5.U7) 


V    >   e 

n,o    °f,o' 


0  <  6  -S  6, 


In  the  region  behind  the  shock  it  follows  directly  from  (5.3)  that 


(^•i^s)   vj^,>bj^,  , 
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0  <  e  <  e, . 


We  shall  now  show  that 

(5.I49)     ^n  1  ^  4,1  '  0  <  e  <  e^^. 

hs   a  first  step  in  this  direction  we  observe  that  (3.2U)  may  be  rewritten  as 
follows  : 

dp,         dH  T       ,2 

where  we  have  made  use  of  the  relation  (3.22)  and  the  fact  that  ^^  n   "  in„/p-.     If 

(5.51)  p^-g(p^,S^) 

is  the  equation  of  state  for  the  medium  behind  the  shock,    then 

dPn  H  .  .       ..  9         Act         ^T 


(5. 


^^)  d^     =     4  ^(PrS,(p,))     .     a,2  .  ^    ^ 


which,  in  virtue  of  (5.ii3)  and  the  fact  that  c  ■  R/(y-1),  becomes: 


V 


"^Pl  2       v-1        clS, 

'5-53)  ^'    a,^  .  X^  P,g^  . 

Thus  we  find  that 

2  dS 

This  equation  may  be  recast  into  a  most  convenient  form  by  making  use  of  the 
following  relations: 


In  the  proof  which  follows  we  shall  mit  the  sub     and  superscript   'f. 
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(5.55) 


(5.56)        iL  = 


sin  e^  V^  T 
0     n^l 

2  2 

n,l         n,l 


(5.^7)    i^  =  ^  ^  w^  * 


RP.      I    ,,o  KHj2   I     ^2 


[^Jil 


dpT       2Pt      1    ^'-i  j^         f     dp,    • 


n 


Equation  (5.55)  is  simply  (3.31)  revrritten  in  a  form  most  suited  to  the  present 
analysis?  equation  (5.5?)  follows  immediately  from  (3.21)  on  making  use  of  the 
ideal  gas  relation  Pn  "^n  =  HT,}  (5.56)  follows  directly  from  the  first  and  last 
expressions  in  (5.3)  after  solving  for  h/^  . 
These  relations  and  (5.5U)  yield 

(5.58)        v|;^j  -  (^2  .  ^,^-h  v2_^  *  b^^,  4  *  (v^^,  -  b^^^,)  ^  "i  =  0. 


where 


(5.59)         G^-  G^(p^) 


1 
TFT 


W' 


L^Hj' 


n 


1-^1 


It  is  important  to  note  that     G,       is  a  positive  quantity  in  fast  magnetic 
shocks  where  [r]  <  0,  [hJ  >  0  and  r^    <     2/(y-1).     Solving  for  V^  ^  in  (5.58) 


we  find  that 
(5.60)     yI  t 

n,l 


_  a^  .  ^.^p,-^  .  /{.I  .  ^.H^p-^)^  -  k^l^4  -  U(V^„^.,  - 


>,2     ^^  ^ 


where  the  choice  of  the  plus  sign  before  the  radical  is  dictated  by  the  requirement 

that  V  ,  approach  c«  _  when  p,  approaches  p  .  Qnploying  the  fact  that 

V  ,  >  b  ,  in  M»  shocks  and  the  fact  that   G,  >  0,  we  find 
n^i    n^i    I  1 
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(5.61) 


2     4-^>4pi'-/4^^pI'^'-^<^ 

n,l  2  '     f?  "*  PQ   • 


It  is  now  a  simple  matter  to  verify  that  the  right-hand  side  of  (5.6l)  is 
c»  ,}  all  that  is  involved  is  a  change  of  notation  [cf.   (l.L)  and  (l,5)j. 


6.     Slow  magnetic  shocks  (h    >  0,  0  <  9     <  90°):     M     shocks 
^  s  o  s  — — — — 

6.1  Slow  magnetic  shocks  of  types  1  and  2;  M^   ^  and  M^    '  shocks 
Our  analysis  of  M     shocks  is  based  on  the  following  relations: 

(6.1)       ^  .  -n  ^  .  -  7 2 a-^  ,  b=    .  /I^   , 

b-V  o  o"-  -J 

n,l  n,l 


(6.3) 


'wS     \2  /„s     \2  2 

n.l  \  1    /     n.o  \  8 


(6.U) 


Hi/         ^^VV;  (P^7;7       [l.(Vh3)sin9j    * 

l+(>i„/h  )sin  »          ^  2  sin  9     +   (r-Dh, 

S        5                      U  OS 

1/h^  -  sin  &                  n^  I (1-s^)  -  X  h  sin  9„  ±  /r  (h^) 

1  -  (Xg/hg)gin  e^           ^s  2  s^  sin  &^  +  (y-Dhg 


(6.6) 


cKTg/hg)  -[(Y-1)(^^A^)2  .  Y  sin  &^  ?3/h^  .  l] 


dh 

s 


i/FoTT 


where 
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d(X  /h  )          2q 
(6,7)       ^     °      o     cos  Qn 


(6.8)  a)  Rjhg)  =  [(1-3^)  -  ^  hgSineJ^  -  (2s^3in  »^+(Y-l)hg)(hg-2  sin  ©^) 

^^  ^x(^s>  -  [i  ^^^"  ®o  -^  (1-^0^1^  -  (2  3in  9^*(Y-l)h3)(h3-2s^3in  9^) 
and 

(6.9)  R_(hg)  -  ^(hg)  =  h^  j  Y-l^nJ^  -  (y-1)1  -  hgsin  e^(2-Y)(l+s^) 

+  (1-3q)  +  USpSin^jj. 

The  simplest  way  to  derive  (6.1)-(6.$)  and  (6.7),  (6.8)  is  to  recognize 

that  (5«l)-(5.5)  and  (5.7),  (5.8)  are  solutions  of  the  shock  relations  for  negative 

as  well  as  for  positive  values  of  h_.  If  we  set 

H   —  H 
(6.1D)    hg  =  -  h^  =  ^^i2l_Z2l  ,        h^.<o, 

o 

and  substitute  (-h  )  into  (5.1)-(5.5),  (5.7),  (5.8)  we  obtain  the  corresponding 
s 

formulas  (6.l)-(6.5),  (6.7),  (6.8).  To -derive  (6.6)  we  substitute  (-h  )  for  h. 

3         X 

in  (5.11),  differentiate  with  respect  to  h  ,  solve  for  ??  /h     and  employ  the 

S  So 

last  equation  of  (6.5)  to  remove  ^_/h  from  the  dencminator  of  the  resulting 

s  s 

expression. 

According  as  the  plus  and  minus  sign  appears  before  the  radical  in  (6.U) 

and  (6,5),  we  shall  speak  of  the  plus  or  minus  branch  of  X  /h  and  ^  /h  and  of 

s     s  s     s 

derived  quantities.     As  before,   the  plus  branches  of  vf  /h     and  X  /h    will  be 

S3  3       3 

denoted  by    "l/h    and  X„A  ,  the  minus  branches  by   v?~A  ,  f/h  .     It  is  easy 
ss  ss  ssss 

to  verify,  using  the  first  equation  in  (6,U),  that  v?  /h     corresponds  to  X  /h 

s     s  s     s 
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and  'L^/hg  'to  ^g/^s*  "^^^   (6,6)  and  (6.7)  we  obtain  the  following  useful  result: 

The  h  derivatives  of  ^   /h  and  X  /h  are  negative  or  positive  according  as  we 

are  dealing  with  the  plus  or  minus  branches;  it  is  assumed  here  that  h  is 

limited  to  those  values  for  which  /R(h  )is   real. 
s 

Employing  a  procedure  similar  to  that  used  in  determining  admissible 
fast  shock  solutions,  we  find  that  admissible  slow  shocks  may  also  be  conveniently 

classified  into  two  types.  Slow  shocks  of  type  1  are  those  in  which 

2  2 

s  >  1  -  Y  sin  9  :  slow  shocks  of  type  2  are  those  in  which  s  <  1-  v  sin  ©  . 
o  o  o  o 

In  the  following  we  shall. give  a  brief  description,  based  on  the  second  equations 
of  (6.1i),  (6,5)  and  the  remark  underlined  above,  of  the  »T  /h  vs.  h  and  X  /h    vs. 

^  S   3        S        S   S 

h  curves  of  types  1  and  2. 
s 

6.2  Description  of  the  curves  ^  /h  versus  h  and  X  /h  versus  h 
's  s  3  s  s  ^— — ^  s 

(l)  2 

6.2.1  M^  '^  curves   (s^  >  1  -  y  sin  9  ) 
S   0  -  o 

Here,  admissible  M  shocks  exist  only  in  the  range 

3 

0  <  h     <  li       =     2  sin  9^   . 

3  3  0 


Within  this  range,  moreover,  only  the  plus  branch  rl  /h     is  admissible.     As 

3   S 

h  increases  from  zero  to  h  ,  ^  /h  decreases  monotonically  [see   (6,6) J  from  the 


s 

finite  value 

s'  *-s     s 

(6,11)       ifVh 

h 

s 

„(l-3^)-H/(l-s^)2*Us^sin2^^ 

"  "                       2s  sin  © 
0            o 

to  zero  [see  (6,5)].   Note  that  H^  ■,   decreases  from  IL.  ^  to  -H 

as  h  increases  from  zero  to  2  sin  ©  .  A  typical  (>7  /h  )  vs.  h  -curve  of 

S  O  3   S         S 

type  1  is  illustrated  by  the  solid  curve  of  Figure  10(a).  The  behavior  of 

X  /h  may  be  ascertained  from  the  behavior  of  ^  /h  using  equations  (6,U)  and 

S3  S   S 

(6.7).  In  particular,  from  (6.7)  it  follows  that  X  /h    varies  in  the  same  sense 

s  s 
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V^/h^  vs    hg 


^(I_s/+4s,sin2«^  +  (l-s„) 
2SoSinfi„ 


-h,  =  2  sin  e„ 
s  o 


2cos«     V(y-i)(i_s  )^+  r^»«S'"^o  -(2-r)('  +  So)sinfi, 


•o'^-'-o 


2,i-2i 


2(y-l)  —  y'=sin^6lo/2 


Xs/hs  vs    h^ 


^(l-s/+4s-sin2e+  ( I  -  s„) 


0'  "0  o 

2  sinfl„ 
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Figures  IDja),   lP(b).  Sketch  of  the  variation  of  ?  /h  xd.th  h  (a), 

S       S  3 

and  of  X  /h    with  h     (b),  in  M     shocks.     The  (±)   superscripts  on  ^     and 
ssss  s 

X     indicate  the  plus  and  minus  branches.     The  solid  line  curves  are  curves 
s 

of  type  1}  the  broken  line  curves  are  of  tjrpe  2.     X  /h     and  ^  /h  , 
whatever  the  type,  vary  in  a  sense  opposite  to  h  j  X~/h     and  ^"/h    vary 

S  3       3  S       S 

in  the  same  sense  as  h  , 

s 
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as  T  Ai  •  The  solid  curve  of  Figure  10(b) illustrates  a  typical  (X  /h  )  vs. 

S   S  S   3 

h  -curve.      In  contrast  to  the  behavior  of  xt/h.,  X*Ai_  is  bounded  for  all  allowable 

S  X       X  S       9 

values  of  h  . 
s 

6.2.2     M^^^  curves     (s^  <  1-^  sin^^) 

M^  shocks  of  type  2  exist  only  in  the  range 


-  sin  eQ(2-Y)(l+s^)  +  2cos  9^/(y-1)(1-Sq)^+  s^Y^sin^©^ 


(6.12)  0  <  h^  <  h^     »     P »r     , 

^  "     ®  2(y-1)  -  (y  3in\/2) 

^  A 

where  h     is  a  zero  of  R(h  ).     In  the  range 
s  s 

(6.13)  0  <  h    ^  h       »     2  sin  e„, 

s         s  o' 

only  the  plus  branch  ^  /h     is  admissible.     In  the  range 

(6.11i)  hg    <    hg    <    hg 

both  plus  and  minus  branches  are  admissible.  On  the  plus  branch  the  values  of 

—  ♦  A 

ri  /h     exceed  those  on  the  minus  branch  except  at  h     =  h    where  the  branches  join 
s     s  s         s 

^      _+ , 

smoothly.  As  h  increases  from  zero  to  h  ,  >?  fh    decreases  monotonically  from 
s  s   s  s 


the  value  ?  A  L  _r>  given  in  (6,11)  to  (>?  A  )  L  ^  •  ^s  h„  decreases  from  h„ 
s  s  h  =0  ^  s  s  h  =h       s  s 

'  s  '  s  s 

to  h  ,  ♦fVh  decreases  monotonically  from  its  maximum  value  (^JA^L  'y,     to  zero, 
3   s  s  s  s  I  n  =n 

The  general  features  of  the  (X  A  )  vs.  h  -curves  are  the  same  as  those  of  the 

3  s      s 

(Jl   A  )  vs.  h  -curves  because  of  (6.7).  /gain  (X  A  )  is  bounded  for  all  allowable 
s  s      s  s  s 

valuesof  h  ,  The  broken  line  cuives  of  Figure  10(a)and  10(b)  illustrate  the  general 

behavior  of  >?  A  and  X  A_  in  the  allowable  range  of  h  . 
s  s     s  s  s 
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6.3  Concluding  remarks 

It  should  be  clear  at  this  p)oint  that  our  analysis  of  M  shocks  may  be 

s 

made  to  parallel  closely  that  of  M„  shocks.  It  is  a  fact  that  almost  every  M„  result 
has  a  parallel  M  result.  Evidently  there  is  no  need  to  give  proofs  of  M  results 
where  this  is  the  case.  Some  M  results,  however,  have  no  parallel  in  M^  results. 
Of  these  the  only  one  which  does  not  follow  more  or  less  directly  from  the  results 
of  the  preceding  two  subsections,  and  which,  we  believe,  merits  special  attention, 
is  the  statement  that  fsl  ,  m  ,  V   ,  when  regarded  as  functions  of  h  ,  reach  their 

S  -At 

maxima  at  the  same  value  of  h  and  ■ti'iat  at  this  point,  V  ,  *=  c  , . 

—————— —^—^——  s  * '  n,l    3,1 

To  prove  this  statement  it  clearly  suffices  to  show  that  V^  ,  =  c  -  at 
^  ''  n,l    s,l 

that  value  of  h  where  m  attains  a  maximum:  for,  m  =  p  V    and,  as  we  have 
s       s  ''so  n,o 

seen  [see  subsection  3*23,  m  varies  in  the  same  sense  as  [sL. 

Now  following  the  procedure  employed  in  subsection  ($.$),  it  is  easy  to 
verify  that 


dm. 


(6.15)   2ivi,f .  (4  ♦  ^^Pih  -/(4  ♦  ^*iPi'^X,i4  -  ""m-'m^^p^  h 

where  the  function  G,   is  given  in  (5.59).  '■.'hen  dm  /dp,  vanishes  we  find,  using 
(l.h),  (1.5),  (6.15)  and  the  fact  that  G^  is  here  Glnite,  that 

(6ol6)     V^  T  =  c^  T 
n,l     s,l 

To  complete  the  proof  it  is  clearly  sufficient  to  prove  l)  that 

dm     dm   dh 
s  ^  s  s 

dp^  "     dhg  dp^ 

vanishes  when  dm  /dh  vanishes  and  2)  that  m  regarded  as  a  function  of  h  has  at 
s   s  s  s 

least  one  stationary  value. 


There  is  only  one  value  of  h  for  which  this  statement  is  truej  but  as  we  have 
mentioned  in  the  footnote  on  p.  22  the  proof  is  too  lengthy  to  be  included  in  this 
report. 
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Statement  1)  will  be  proved  if  it  can  be  shown  that  dm  /dh  and  dp^/dh 
do  not  vanish  simultaneously.  Now  it  is  easily  verified  that 

(6  17)   A.  ^^"s/Po^  ^  [l-  sin  »,/h^](d?3/dh^)  *  \   (^3^l)sin  9^Ag 
hhJ    ^s    "  [l*(?/h^)sinej2 

Since  ©^  t'  0,  h„  >  0  and  7^  >  0  it  is  clear  that  dm„/dh„  and  d'i?„/dh„  cannot 
o''s        's~  ss      ss 

vanish  simultaneously. 

Statement  2)  is  an  immediate  consequence  of  Rolle's  Theorem.  Consider 

M    shocks  for  example.  In  this  type  of  shock 
s 


(6.18)   /  = 


-1 1/2 


n,o 


l+(?3/h^)sin»^ 


(•'^s-^^^c 


l*{l[^\)si.n^^ 


1/2 


[cf.   (6.3)]   is  a  single-valued  differentiable  function  of  h  ,  0  <  h    <  2  sin  d  , 

which  assumes  the  value  b        when  h     »  sin  ©     and  when  h     ■  2  sin  &„  (Note  <?     "0 

n,o  s  o  s  o  *3 

at  h    =  2  sin  ©  ):  dm  /dh    therefore  vanishes  at  some  interior  point  of  the 
s  o  '       s       s 

(2) 
interval  sin  ©     <  h    <  2  sin  ©  .     Similar  reasoning  applied  to  M        shocks  leads 

0-S~  O  Ofl-  g 

to  the  conclusion  that  dm  /dh  vanishes  for  at  least  one  interior  point  of  the 

s   s 

interval,  sin  8  <  h  <  h  ,  where  h  is  given  in  (6.12),  A  more  detailed 
o  ~  s  ~  s        s 

analysis  reveals  that  the  vanishing  dm  /dh  occurs  at  a  value  of  h  exceeding 

s   s  s 

1.?  sin  ©„. 
o 


This  equation  follows  from  (6.3)  and  the  identity 


d7  /h 
s'  s 

~aE — 


1   ^^S    ^^3 
S 
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7.  Limit  shocks 

7.1     Fast  0°-limit  shocks 

7.1.1  Fast  0°  -  limit  shocks  of  type  1  (s^  ^  1) 

The  analysis  of  this  subsection  is  based  on  the  identity 

d"^f    7f       d(7-/h  ) 


and  the  equations  [cf.  (5.5)  and  (5. 8) J. 

2  s^sin  e^  -  (Y-l)h, 


(7.2) 


^f   Y 


-^  h^sin  ^^   -  (1-3^)+  /^(hj) 


(7.2«)    R(h 


where 
(7.3) 


If)  =  [j-j       \   Y^3in^^-U(Y-1)  ■  ♦  hfSin  &^{2'^){l*s^) 


+  (1-s^)^  +  h   s^sin^©^. 


2  sin  e. 


^    <   ^f   ^   "yIT—  ' 


and  where  only  the  plus  sign  appears  before  the  radical  in  accordance  with  the 
fact  that  we  are  dealing  with  shocks  of  type  1, 

Now  let  us  suppose  that  n ,  0  <  rj  <  2/(y-1),  is  fixed  and  let 
hf  =  h„(^-,&  )  bfe  the  corresponding  value  of  h_.  It  follows  immediately 
from  (7«3)  that 


(7.U) 


» 


■^0  VV®o^  "  °'   °<V  Y^ 


Lim 


If  we  take  the  limit  of  both  sides  of  (7.2),  we  find  that  non-vanishing  values  of 
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7  ,0  <n^<   2/(y-1)j  are  compatible  with  the  vanishing  of  h^  in  the  9^  =  0 
limit.  Employing  (7.1),  (5.6)  and  (5.7)  we  find 

dh- 
(7.5)     Lim   — i   =0,     0     <      K^  <     2/(y-1). 

Thus,  in  the  (h_,  ?«)  plane  the  limiting  process  may  be  described  as  follows: 
Ps  9     approaches  zero  the  range  of  h,.  shrinks  and  the  tf^  vs.  h^- curves  become 
steeper.     In  the  limit  h-  vanishes  and  ^-  becomes  independent  of  h^. 

Taking  the  limit  of  both  sides  of  equations  (U.5-B^-B,  )    we 
find,  in  virtue  of  (7.U),  that 


2?fir 

2-(y-1)9. 


(7.6)  Y^  "     \   /  _^_       ,       ^<'if<  2/(y-1)  > 


(7.7)  m^     =     -[p]f/OTf, 

(7.8)  Wf  =  -'^f'^^f» 

(7.9)  [uylf     =    0. 

These  equations  will  be  recognized  as  the  well-known  gas  dynamical  Jump 
relations!  they  lead  directly  to  the  fast  gas  shock  equations,  namely  (2.U5). 

7.1.2     Fast  0°-liinit  shocks  of  type  2  (s^  <  l) 

Consider  the  A.^  vs.  h-  curves  of  Figure  3(tf'  Note  that 
h^  =  h„('T-,  9   )  approaches   zero  as  9     approaches  0     on  the  upper  branch  (minus 
branch)  in  scwie  neighborhood  of  ^  =  2/(y-1)  while  outside  this  neighborhood  h^ 
approaches  a  finite  value  depending  on  rf  „   [see  0  -curve  of  Figure  3(ty] . 
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These  observations  may  be  verified  analytically  as  follows.     Let  us  take  the 
limit  of  both  sides  of 


(7.10)    ^  = 


^f      [-  i  V^  ^0  -  (^-^o^  -  /^^^  ] 


f        2sQsin  B^-   (y-1)  h^ 


where 


(7.n)   0  <  h^  <  ^^(e^)  , 

[see  (5.5)  and  (5.23)J  without  assuming  that  h„  also  approaches  zero.  We  then 
find  that 

(1-s  )  1  /(1-s  )^  -  (Y-l)h? 


(7.12)   n^ 


(Y-1) 


where 

(7.13)   0  <  h„  <  (1-s.)/  /Tr^. 


2 

Equation  (7.12), when  solved  for  h„  ,  becomes 


^f '  -  7f  [ 


(7.11:)    h/     "    ff^  f  2(l-s^)  -  (Y-l)*Jf 


2 

The  condition  h-  >  0  implies 

2(l-s  ) 

(7.15)   0  <  ^^  <  -^:^   =  ^^^^^.^^  <  2/(Y-i). 

Assuming  that  ^„  is  confined  to  this  range,  and  tdcing  the  limits  of  both  sides  of 
(5.1),  (5.2),  (5.3),  (5'.U)  we  obtain  a  set  of  equations  which  lead  directly  to 
those  listed  in  (2.5.1)  -  that  is,  to  the  equations  which  determine  the  state 
behind  an  ' incorplete '  switch-on  shock.   We  now  ask  what  is  the  limit  solution 
corresponding  to  values  of  7^     in  the  range 
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(7.16)  ^f,cri.t^'^f-  2/(y-1). 

To  determine  this  solution  we  expand  (7.10)  in  the  neighborhood  of 

h-  =  2  sin  ©  /(y-1)  and  assume  that  h„  is  small  when  &  is  small.  The 

result  is 

2(1-5  ) 

(7.17)  nr 


(y-1)-2  s^sin  e^/h^ 

where  the  dots  indicate  terms  of  higher  order  in  h-  and  sin  9  .  From  (7.17) 
it  follows  that 

2(1-3  ) 

(7.18)  Pf,  >  -^^  =if^^^^,, 

for  any  small  value  of  h^  and  sin  &  ,  Solving  for  h»,  we  find 

23.sin  9 

(7.19)  h-  =  —2 2 ;f  +  ...  . 

'     (Y-l)(?r  ^,,erit^ 
From  this  equation  it  follows  that 

(7.20)  Limhf=  Lijn  ^^(yi^,^^)     "    0>    ^f^crit  "^  ^f  -  7^  * 
o"      o" 

By  following  the  procedure  employed  in  subsection  (7.l)>  it  is  now  an  easy 
matter  to  verify  that  the  limit  solution  is  an  'incomplete'  fast  gas  shock 
[see  Subsection  2,6,2].  It  is  also  a  simple  matter  to  verify  that  the  state 
behind  the  'incomplete'  switch-on  shock  passes  continuously  into  the  state 
behind  the  incomplete  gas  shock  as  •?_  passes  through  the  critical  value,  ?gj^4.« 


-  91  - 


7.2  Fast  90°-llmit  shocks  (s  >  l) 

We  have  seen  that  90  -limit  shocks  are  all  of  type  1.  If  we  take  the 

limit  as  9  approaches  90  of  both  sides  of  the  last  equations  in  (5.U)  and 

(5.5)  we  find,  using  the  fact  that  only  the  plus  branch  is  admissible  in  type  1 

shocks  and  the  fact  that  s  >  1, 

0  - 


(7.21)  ri^/h^    =  1. 

This  equation  together  with  the  limit  equations  resulting  from  (5.1)  -  (5.3) 
and  the  Hugoniot  relation  lead  directly  to  the  equations  listed  in  (2.65)  - 
that  is,  to  the  equations  which  determine  the  state  behind  fast  perpendicular 

shocks. 

7.3  slow  limit  shocks 

The  various  slow  limit  shocks  described  in  subsections  2.6.3  and  2,6,5 
may  be  obtained  as  limits  of  (6.1)  -  (6.5)  or  (U.5)  B^  -  B,  by  employing 
methods  similar  to  those  used  abovej  a  detailed  derivation  will  therefore  be 
omitted. 
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